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Abstract 


We  use  mathematical  modeling  in  decision  analysis  to  help 
us  obtain  a "better”  profit  lottery  than  we  can  assess  directly. 

The  concept  of  the  authenticity  of  probabilities  is  introduced  to 
define  the  measure  of  "goodness"  of  the  profit  lottery.  The  role 
of  modeling  is  to  simplify  our  assessment  task  through  the 
decomposition  of  the  profit  lottery.  However,  budgetary  constraints 
force  us  to  make  approximations  in  the  modeling  process  and  thereby 
cause  us  to  misstate  the  profit  lottery.  The  models  used  in  a 
decision  analysis  should  be  regarded  as  subjective  expressions  of 
our  uncertainty  rather  chan  as  objective  descriptions  of  the 
real-world. 

A methodology  is  presented  that  quantitatively  relates  the 


modeling  approximations  made  in  a decision  analysis  to  the  results 
of  the  analysis. 
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INTRODUCTION 


To  the  Wizard  and  to  the  Clairvoyant,  decision  making  is  a 
simple  matter.  The  omnipotent  Wizard  has  the  power  to  change 
anything  in  the  world  to  suit  his  own  taste.  To  him,  there  is  no 
such  thing  as  a decision  problem;  he  merely  causes  to  happen  what- 
ever he  wishes.  The  omniscient  Clairvoyant,  although  lacking  the 
Wizard’s  power  to  change  events,  possesses  the  power  to  foretell 
the  future  perfectly.  To  him,  there  is  no  uncertainty  in  decision 
®^king , he  merely  chooses  the  course  of  action  whose  consequences 
he  most  desires. 

By  contrast,  we  lowly  mortals  cannot  pretend  to  possess  the 
Wizard's  or  the  Clairvoyant's  powers  in  decision  making.  But 
there  is  a superbeing  whom  we  do  try  to  emulate  — the  Elicitor. 
The  Elicitor  can  neither  change  nor  foretell  events  and  he  is 
therefore  uncertain  about  the  future.  However,  he  possesses  the 
special  ability  to  fully  and  accurately  express  his  uncertainty  in 
the  form  of  probability  statements.  In  a decision  problem,  the 
Elicitor  assesses  his  probabilities  on  the  future  consequences  of 
the  alternative  courses  of  action  and,  using  his  utility  function, 
calculates  the  expected  utility  for  each  alternative.  He  then 
chooses  the  course  of  action  having  the  highest  expected  utility. 

He  is  not  assured  that  his  decision  will  lead  to  the  most  desired 
outcome,  but  he  is  confident  that  his  actions  are  wholly 
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consistent  with  his  preferences  and  with  his  uncertain  under- 
standing of  the  future. 

In  decision  analysis,  seeking  to  imitate  the  Elicitor,  we 
likewise  make  probability  statements  about  the  future  consequences 
of  our  actions  and  calculate  the  resulting  expected  utility  for 
each  alternative.  But  we  unfortunately  do  not  possess  the 
Elicitor’ s ability  to  directly  assess  probabilities  that  fully 
and  accurately  express  our  uncertainty.  Recognizing  this,  we  use 
modeling  to  help  us  obtain  the  requisite  probability  statements. 

Modeling  is  the  source  of  most  of  our  dissatisfaction  with 
the  results  of  particular  decision  analyses.  While  we  do  not  con- 
test the  validity  of  decision  theory,  we  often  complain  that  the 
models  used  in  a decision  analysis  are  "too  simplistic"  or  "not 
realistic  enough"  or  "not  believable"  and  we  therefore  regard  the 
results  of  the  decision  analysis  with  doubting  eyes. 

How  should  we  deal  with  our  dissatisfaction  about  modeling  in 
decision  analysis?  What  do  we  mean  by  "goodness"  in  a model  and 
can  we  quantify  it?  Can  we,  for  instance,  define  an  index  of 
"realism"  or  of  "credibility"  on  models?  How  should  we  choose 
among  alternative  models?  For  example,  can  we  use  the  notion  of  a 
probability— space  of  models?  And  how  should  we  decide  when  to  do 
more  modeling?  For  instance,  could  we  use  the  concept  of  the  value 
of  perfect  modeling  as  an  analogue  to  the  value  of  perfect 
information? 
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In  this  dissertation,  I offer  a way  to  think  about  modeling 
in  decision  analysis  so  that  we  can  deal  with  our  dissatisfactions 
meaningfully.  Chapter  1 provides  a philosophical  perspective  on 
the  role  of  modeling  in  decision  analysis.  Chapter  2 presents  a 
methodology  with  which  we  can  quantify  our  dissatisfaction  about 
the  modeling  in  a decision  analysis  and  relate  it  to  the  results 
of  the  analysis.  Chapter  3 is  an  example  illustrating  the  use  of 
the  methodology  and  Chapter  4 is  an  extension  of  the  methodology 
to  stochastic  models. 
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CHAPTER  1 

A PERSPECTIVE  ON  MODELING  IN  DECISION  ANALYSIS 


1 • 1 Introduc  t ion 

In  theory,  we  do  not  need  to  use  mathematical  modeling  to 
perform  decision  analysis.  Analytically,  a decision  problem  is 
defined  by  two  variables:  a decision  variable  d , represent- 

ing the  alternative  courses  of  action  open  to  us,  and  an  outcome 
variable  v , representing  the  relevant  consequences  of  our 
actions.  In  decision  analysis  [7,8]  , we  need  only  assess  for 

each  value  of  d the  conditional  probability  distribution  on  v 
given  d , {v|d,4}  , called  the  profit  lottery  and  state  our 

risk  preference  oy  specifying  a utility  function  on  v , u(v). 

Then,  we  can  calculate  the  expected  utility  of  each  profit 
lottery: 

d lv  {v|d,4}  u(v) 

and  identify  that  value  d*  of  the  decision  variable  correspond- 
ing to  the  profit  lottery  with  the  highest  expected  utility  (i.e., 

the  most  preferred  profit  lottery).  The  optimal  alternative  is 
the  one  specified  by  d*. 

Although  superfluous  in  theory,  modeling  is  nevertheless 
indispensible  to  decision  analysis  in  practice.  In  theory,  we 
assume  that  we  can  directly  assess  the  profit  lottery,  but  in 
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practice,  we  are  unable  to  do  so  satisfactorily.  Consequently, 
we  turn  to  modeling  to  help  us  obtain  the  profit  lottery. 

Our  ideal  in  decision  analysis  is  not  to  construct  the 
perfect  model,  but  rather  to  obtain  the  authentic  profit  lottery  — 
the  one  that  accurately  expresses  our  uncertainty  about  the  future. 
(See  Figure  1.1  for  a conceptual  "roadmap".)  We  do  not  have  the 
Elicitor's  ability  to  assess  the  authentic  profit  lottery  directly, 
so  we  employ  the  modeling  strategy,  which  would  yield  the  authentic 
profit  lottery  if  successfully  executed.  However,  budgetary  con- 
straints force  us  to  make  approximations  in  the  modeling  process ; 
consequently,  the  profit  lottery  that  we  obtain  through  modeling  is 
not  the  authentic  profit  lottery.  The  modeling  approximations, 
then,  are  the  sole  source  of  our  dissatisfaction. 

1.2  Authentic  Probabilities 

We  use  modeling  in  decision  analysis  because  we  do  not  believe 
that  the  directly  assessed  profit  lottery  is  good  enough. 

But  what  is  "goodness"  in  a profit  lottery?  The  profit  lottery, 
as  a subjective  probability  statement  |j3,5,16^  , is  the  quantified 

expression  of  our  beliefs  about  the  likelihood  of  occurrence  of  real- 
world  events  beyond  our  immediate  perception.  A good  probability 
statement  is  simply  one  that  accurately  and  fully  expresses  our 
beliefs.  To  denote  a good  probability  statement,  I use  the  term 
"authentic",  which  the  dictionary  defines  as  'Worthy  of  acceptance 


and  belief." 


Authentic 


limitations 


profit 


Budgetary 


Modeling 


limitations 


Modeled 

profit 

lottery 


Modeling 


Authenticity  is  not  the 


snme  as  truenebs.  The  only  test  of 
goodness  of  a probability  statement  is  whether  or  not  we  believe 
it,  not  whether  it  is  true  or  false.  Indeed,  it  is  not  even 

meaningful  to  speak  of  a probability  statement  as  being  true  or 
false. 

Compare,  for  example,  the  two  statements:  "It  will  rain 

tomorrow.  and  The  probability  of  rain  tomorrow  is  3/4."  The 
first  statement  is  an  assertion  about  the  real-world  and  is  either 
true  or  false.  We  can  determine  its  goodness  (i.e.,  its  trueness) 
by  observing  the  world  (tomorrow).  The  second  statement,  on  the 
other  hand,  does  not  describe  the  real-world,  but  instead 
expresses  our  uncertainty  about  the  world.  It  makes  no  sense  to 

83y  that  lt:  is  true  or  false,  but  only  that  it  is  authentic  or 
inauthentic.  We  cannot  ascertain  its  goodness  (i.e.,  its  authenti- 
city) by  appealing  to  the  real-world,  but  rather  only  to  our  beliefs, 
The  authenticity  of  a probability  statement  depends  on  our 
state  of  information.  A probability  statement  that  we  accept  as 
authentic  based  on  one  set  of  information  may  become  inauthentic 
when  we  receive  additional  information.  For  example,  if  „e  observe 
rainclouds  gathering  overhead,  our  authentic  probability  for  rain 
tomorrow  may  change. 

But  how  do  we  obtain  authentic  probabilities?  Clearly,  the 
authentic  probability  for  «ja  event  must  reflect  all  of  the  informa- 
tion we  possess  relevant  to  that  event.  Therefore,  to  assess  the 
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authentic  probability,  we  must  review  our  entire  state  of 
information  pertaining  to  the  event  and  judge  the  likelihood  of 
its  occurrence  on  the  basis  of  that  information. 

This  is  a formidable  task  for  almost  any  event.  Since  we 
believe  that  the  world  is  highly  interdependent,  our  state  of 
information  ic  an  extensive  and  complicated  web  of  knowledge. 
Hence,  we  must  retrieve,  organize  and  process  a large  amount  of 
information  to  assess  an  authentic  probability. 

Consider,  for  example,  the  probability  that  the  winner  of 
the  1980  U.S.  Presidential  election  is  a member  of  the  Republican 
Party.  To  assess  the  authentic  probability  for  this  event,  we 
must  consider  every  set  of  circumstances  that  might  lead  to  its 
occurrence.  We  must  review  our  knowledge  of  historical  trends, 
of  the  current  state  of  the  Republican  Party  and  of  the  intentions 
and  qualifications  of  potential  candidates.  We  must  relate  the 
authentic  probability  to  our  uncertainty  about  intervening  events, 
such  as  the  outcome  o'  the  197C  election  and  the  possibility  of 
war  or  depression  before  1980. 

We  encounter  two  major  difficulties  in  trying  to  assess  the 
authentic  probability.  First,  we  find  that  the  task  of  consider- 
ing everything  that  might  affect  the  outcome  of  the  event  is 
seemingly  endless.  And  second,  we  find  it  virtually  impossible 
to  perform  the  requisite  processing  of  the  information  without 
external  computational  assistance. 
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Let  us  Imagine  the  existence  of  a person,  called  the 
Probabilist,  who  is  capable  of  performing  upon  request  any 
calculation  using  the  rules  of  probability  calculus  (e.g.,  Bayes’ 
Theorem,  expansion,  change  of  variable1* . Then,  we  can  state  the 
following  operational  definition  of  authenticity:  The  authentic 

probability  for  en  event  is  the  one  we  would  obtain  if  we  could 
spend  an  unlimited  amount  of  time  in  introspection  and  if  we  had 
available  the  services  of  the  Probabilist. 

It  should  be  readily  apparent  that  the  authentic  probability 
is,  in  most  cases,  an  unattainable  ideal.  Even  after  lengthy 
introspection,  we  can  almost  always  think  of  something  relevant  that 
we  have  not  yet  considered.  And,  of  course,  we  do  not  possess  the 
computational  capabilities  of  the  Probabilist. 

1.3  Operative  Probabilities 

Since  we  cannot  obtain  authentic  probabilities,  we  must  use  the 
probabilities  that  are  based  on  only  partial  consideration  of  our 
state  of  information.  These  probabilities  I call  "operative". 

As  a matter  of  notation,  we  can  represent  our  Incomplete 
consideration  of  the  state  of  information  & as  C(4).  Then,  our 
authentic  and  operative  probabilities  for  event  E are: 

{E | -& } ■ p Authentic  probability 

{e|C(4)}  - p^  Operative  probability 
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It  is  useful  here  to  distinguish  between  two  kinds  of 
uncertainty.  The  first  kind,  called  primary  uncertainty,  is  the 
uncertainty  due  solely  to  the  finiten^ss  of  our  state  of  informa- 
tion and  is  expressed  by  authentic  probabilities.  The  second 
kind,  called  secondary  uncertainty,  is  the  additional  uncertainty 
due  to  the  incomplete  consideration  of  our  state  of  information. 
Whereas  primary  uncertainty  can  be  resolved  only  by  receiving  more 
information,  secondary  uncertainty  can  be  resolved  by  further 
introspection  and  calculation. 

Conceptually,  we  can  represent  our  secondary  uncertainty  as 
a conditional  probability  distribution  on  the  authentic  probability 
yiven  partial  consideration  of  our  state  of  information: 

{pJCOS)} 

a 

Then,  the  mean  of  this  distribution  is  the  operative  probability: 

P0  =<pJc(4)> 

and  the  variance  v(p  |c(4))  is  a measure  of  how  much  secondary 
uncertainty  remains  to  be  resolved  and  thus  how  "close"  the  opera- 
tive probability  is  to  the  authentic  probability. 

As  we  spend  more  time  i:  introspection  and  calculation,  we 
reduce  the  amount  of  secondary  uncertainty  and  thereby  make  the 
operative  probability  converge  to  the  authentic  probability. 
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For  example,  the  operative  probability  of  a Republican 
victory  in  the  1980  election  that  we  assign  after  only  a moment’s 
thought  is  ‘'far"  from  the  authentic  probability;  that  is,  since 
much  secondary  uncertainty  remains,  the  operative  probability  may 
change  significantly  with  further  introspection  and  calculation. 

On  the  other  hand,  the  operative  probability  that  we  assign  to 
getting  a head  on  the  next  flip  of  a coin  is  very  "close"  to  the 
authentic  probability  because  little  secondary  uncertainty  remains 
and  further  consideration  of  our  state  of  information  is  unlikely 
to  markedly  change  the  operative  probability;  in  this  case,  we 
can  say  that  the  operative  probability  is  virtually  authentic. 

But  what  difference  does  it  make  how  "close"  the  operative 
probability  is  to  the  authentic  probability?  If  we  must  act  with- 
out further  consideration  of  our  state  of  information,  then  it 
makes  no  difference  how  "close"  the  operative  probability  is  to  the 
authentic  probability.  However,  if  we  are  able  to  further  consider 
our  state  of  information  before  acting,  then  the  "farther"  the 
operative  probability  is  from  the  authentic  probability,  the 
greater  the  value  of  doing  some  more  introspection  and  calculation. 

For  example,  suppose  that  we  are  offered  a choice  between  a 
sure  $60  and  a lottery  yielding  $100  if  a specified  event  occurs 
and  nothing  otherwise;  the  event  may  be  either  a Republican  victory 
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in  the  1980  election  or  a head  on  the  next  flip  of  a coin. 
Suppoca  further  that  our  operative  probability  pQ  of  a 1980 
Republican  victer>  is  identical  to  our  operative  probability  q 

c 


of  a head  on  the  next  flip  of  the  coin,  but  that  pQ  is  "far* 


from  the  authentic  probability  while  q is  virtually  authentic: 


{plC<4)} 

a. 


1 0 < p < 1 

— a — 


otherwise 


{qa|c<*)} 


6(qa  - >s) 


Then,  if  we  are  prohibited  from  doing  further  introspection  or 
calculation,  it  does  not  matter  which  event  is  used  in  the  Lottery; 
in  either  case,  the  expected  value  of  the  lottery  is  $50  and  our 
decision  is  to  choose  the  sure  $60.  However,  if  we  are  permitted 
the  opportunity  to  introspect  and  calculate  before  choosing,  it  does 
matter  which  event  is  used  in  the  lottery.  If  the  lottery  depends 
on  the  coin  flip,  no  amount  of  introspection  or  calculation  vill 
change  the  expected  value  of  the  lottery,  so  we  would  always  take 
the  sure  $60.  In  this  case,  further  introspection  and  calculation 
has  no  value.  On  the  other  hand,  it  does  have  value  if  the  lottery 
depends  on  the  1980  election.  In  this  case,  further  introspection 
and  calculation  may  reveal  that  the  authentic  probability  is  such 
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that  we  would  choose  the  lottery  instead  of  the  sure  $60 
(Pfl  >»6).  The  expo -ted  value  of  the  decision  given  complete 
introspection  and  calculation  ia: 

v * | <v|Pa>{Pa|C(4)}  dpa 

^a 

- j ’6  6°  dpa  + f 1 100  pa  dp 

o *6 

“36  +32 

- $68 


Therefore,  the  expected  value  of  complete  introspection  and 

calculation  (i.e.,  the  value  of  obtainlt  !>  the  authentic  probability) 
is  $68  - $60  - $8. 

The  value  of  further  consideration  of  our  stace  of  information 
is  analogous  to  the  value  of  receiving  additional  information.  In 
both  cases,  the  value  depends  on  how  much  uncertainty  remains  to  be 
resolved  and  on  how  much  is  at  stake  in  the  decision. 


1.4  The  Modeling  Strategy 

In  decision  analysis,  then,  we  want  the  authentic  profit 
lottery  — the  one  and  only  probability  distribution  that  accurately 
expresses  our  beliefs  about  the  future  consequences  of  our  actions. 
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But  in  most  decision  situations,  otir  state  of  information  about 
the  decision  environment  is  much  too  complicated  for  our  simple 
human  minds  to  consider  completely. 

For  example,  suppose  that  our  decision  is  whether  or  not  to 
manufacture  and  market  a new  product.  We  might  be  able  to  identify 
scores  of  factors  that  affect  the  outcome  of  our  decision,  including 
the  behavior  of  our  competitors,  the  general  economic  conditions, 
the  effectiveness  of  our  advertising,  the  costs  of  raw  materials, 
the  efficiency  of  the  manufacturing  process  and  the  regulatory 
behavior  of  the  government.  The  authentic  profit  lottery  must 
express  our  uncertain  understanding  of  how  each  of  these  factors 
affects  our  profit  and  how  each  will  behave  in  the  future. 

Clearly,  we  would  be  hard  pressed  to  process  all  of  this 
information  mentally.  Consequently,  the  profit  lottery  that  we 
assess  directly  is  very  "far"  from  the  authentic  profit  lottery  and 
the  value  of  resolving  our  secondary  uncertainty  is  high. 

Recognizing  this , we  turn  to  modeling  to  help  us  resolve  our 
secondary  uncertainty.  In  essence,  the  modeling  strategy  is  one  of 
divide-and-conquer ; in  modeling,  we  decompose  the  profit  lottery 
into  smaller  pieces  that  our  simple  minds  can  handle  comfortably. 
Then,  we  rely  on  external  means  of  calculation  to  reassemble  the 
pieces  to  obtain  the  profit  lottery. 
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The  modeling  strategy  is  as  follows: 

1*  Identify  a set  of  real-world  factors  on  which  profit 
is  believed  to  depend,  representing  them  as  state  variables, 
denoted  by  the  vector  s. 

2.  Encode  our  uncertain  understanding  of  the  dependence 
relationship  between  profit  v and  the  decision  and  state 

variables  d and  s as  a conditional  probability  distribution 
{v|d,s,4}  . 

3.  Encodi  our  uncertainty  about  the  future  behavior  of  the 
state  variables  as  a probability  distribution  {s|d}  . 

4.  Using  external  means  of  calculation,  determine  the 
profit  lottery  via  the  expansion  equation: 

{v|d,4}  = [ {v|d,s,d}{s|4} 

_ 

Note  that  the  modeling  strategy  does  not  relieve  us  of  the 
necessity  to  assess  probabilities.  Rather,  it  merely  substitutes 
for  the  direct  assessment  of  the  profit  lottery  the  assessment  of 

probabilities  on  the  state  variables  and  on  the  dependence  relation- 
ship. 

The  key  to  success  of  the  modeling  strategy,  then,  is  select- 
ing a set  of  state  variables  such  that  we  can  handle  the  resulting 
assessment  task  satisfactorily. 
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If  we  can  assess  virtually  authentic  probabilities  on  the 
state  variables  and  on  the  dependence  relationship,  then  the 
modeling  strategy  would  yield  the  authentic  profit  lottery, 
since  the  expansion  equation  is  tautologically  true. 

1.5  The  Abridged  Modeling  Strategy 

Unfortunately,  we  can  seldom  fully  execute  the  modeling 
strategy  because  we  are  limited  by  economic  constraints  on  the 
decision  analysis  Itself.  If  a decision  analysis  is  to  be  worth- 
while, its  cost  must  be  small  relative  to  the  resources  allocated 
in  the  decision.  Therefore,  in  any  decision  analysis,  we  must 
strictly  limit  the  tine  and  effort  devoted  to  obtaining  the 
profit  lottery. 

Full  execution  of  the  modeling  strategy,  by  contrast,  would 
require  much  time  and  effort.  Generally,  the  assessment  of 
virtually  authentic  probability  distributions  on  the  state  variables 
and  on  the  dependence  relationship  would  be  very  time-consuming. 
Furthermore,  because  we  cannot  analytically  perform  integration 
over  an  arbitrary  continuous  function,  as  required  by  the  strategy, 
we  would  need  to  perform  the  integration  numerically,  which  would 
also  be  exceedingly  expensive. 
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Recognizing  that  we  are  unable  to  fully  execute  the  model- 
ing strategy  within  the  budgetary  constraints,  we  abbreviate  it. 
The  abridged  modeling  strategy  is  as  follows: 

1*  Identify  a set  of  real-world  factors  on  which 
profit  is  believed  to  depend,  representing  them  as  state 
variables,  denoted  by  the  vector  s. 

2.  Identify  a deterministic  function  g that 
approximates  the  dependence  relationship  between  profit 
v and  the  decision  and  state  variables  d and  s : 

v a g(d,sj) 

The  function  g is  commonly  called  the  "model", 

3.  Approximately  encode  our  uncertainty  about  the 
future  behavior  of  the  state  variables  as  a discrete 
probability  function  p(£>. 

4.  Using  external  means  of  computation,  calculate 
the  profit  lottery  via  the  expansion  equation: 

(v | d,4 ^ = l 6(v  - gCd.fjO)  p(s) 

£ 

The  abbreviation  of  the  modeling  strategy  results  in  the 
familiar  decision  tree  representation  of  the  decision  problem. 
[23^  (See  Figure  1.2)  The  discrete  steps  of  the  state  variables 
are  represented  as  "branches"  at  the  chance  nodes  of  the  tree  and 
the  profit  assigned  to  each  "tip"  of  the  tree  is  given  by  the 
deterministic  function: 

v = g(d,s) 
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The  degree  to  which  we  must  restrict  the  number  of  state 

variables  and  the  level  of  discretization  of  their  distributions 

depends,  of  course,  on  the  budgetary  constraints  of  the  particular 

decision  analysis.  Usually,  the  restriction  la  quite  severe 

the  size  of  the  decision  tree  grows  geometrically  with  the 

number  of  variables  and  with  the  number  of  discrete  steps  In  the 

distributions,  .tore  specifically,  If  we  have  N state  variables 

and  If  we  discretize  the  i-th  state  variable  into  m±  steps,  the 

number  of  tips  k In  the  decision  tree  for  each  decision  alterna- 
tive  is  given  by: 


k = 


N 

n m. 

i-1  1 


For  each  of  the  k tips,  we  must  evaluate  the  function  v - g(d,s) 

to  determine  the  profit  associated  with  that  tip  and  we  must  multiply 

together  N probabilities  to  find  the  probability  corresponding  to 
that  value  of  profit. 

Clearly,  the  amount  of  computation  required  to  solve  a decision 
tree  can  easily  become  prohibitively  large.  For  example.  If  we  have 
ten  state  variables  and  if  each  is  discretized  into  five  steps,  our 
decision  tree  will  have  nearly  ten  million  tips  for  each  decision 
alternative.  Solving  this  tree  would  tax  the  capabilities  of  almost 
any  computer  and  would  Incur  an  exceedingly  high  cost. 
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Consequently,  to  meet  the  budgetary  constraints,  we 
deliberately  ignore  our  uncertainty  about  some  of  the  state 
variables  in  the  interest  of  computational  economy.  We  select 
from  the  state  variables  only  those  very  few  that  have  the 
greatest  effect  on  profit,  as  revealed  by  sensitivity  analysis, 
and  designate  them  as  aleatory  variables.  Ql9 ,27^]  We  assess  the 
probability  distribution  on  each  aleatory  variable  and  discretize 
it  into  only  a few  steps.  We  designate  the  remaining  state 
variables  as  non-aleatory  and,  disregarding  our  uncertainty  about 
them,  fix  each  of  them  at  a single  representative  value.  In  effect, 
we  discretize  each  non-aleatory  variable  into  just  one  step  instead 
of  into  several.  Typically,  we  have  fewer  than  ten  aleatory 
variables  and  we  rarely  discretize  a distribution  into  more  than 
three  steps. 

The  profit  lottery  that  we  obtain  through  the  abridged  modeling 
strategy  is  not  the  authentic  profit  lottery.  The  profit  lottery 
from  the  abridged  modeling  strategy  isi 

{v|d,4}m  - l 6(v  - g(d ,_s))  p(s) 

£J_ 

From  the  full  modeling  strategy,  the  authentic  profit  lottery  is: 

{v|d,4}  - f {v|d,£,4>  {s_|-&} 

-’s 


20 


The  discrepancy  between  the  modeled  and  authentic  profit  lotteries 
is  due  solely  to  the  two  approximations  made  in  the  abridged 
modeling  strategy: 

1.  {v|d,£,4}  = 6(v  - g(d,s_)) 

2.  {sJ<S}  - p(s) 

I call  the  first  the  Dependence  Approximation  and  the  second  the 
Distribution  Approximation.  In  Chapter  2,  I present  a methodology 
for  estimating  the  size  of  the  discrepancy  in  the  profit  lotteries 
caused  by  these  two  modeling  approximations. 

1.6  The  Subjective  Nature  of  the  Model 

In  the  abridged  modeling  strategy,  we  attribute  all  of  our 
uncertainty  about  profit  to  our  uncertainty  about  the  state  variables 
and  none  to  uncertainty  about  the  dependence  relationship  between 
profit  and  the  state  variables. 

This  arrangement,  although  motivated  by  budgetary  considerations, 
is  also  conceptually  appealing  because  it  mirrors  the  way  in  which  we 
visualize  the  world.  According  to  our  Western  world-view,  we  believe 
that  phenomena  do  not  occur  spontaneously,  but  are  caused  or  influ- 
enced by  other  phenomena;  hence,  we  "explain"  the  occurrence  of  a 
phenomenon  by  establishing  the  occurrence  of  its  antecedents.  Further- 
more, we  believe  that  the  interactions  among  the  phenomena  are 
goverened  by  permanent  rules  of  relationship  (i.e.,  natural  "laws"). 


It  Is  the  role  of  science  to  Increase  our  understanding  of  the 
world  by  discovering  and  codifying  these  rules  of  relationship 
through  the  careful  observation  of  the  universe. 


Applying  this  world-view  to  the  decision  situation,  we 
believe  that  a particular  level  cf  profit  does  not  simply  happen, 
but  is  caused  by  certain  extraneous  factors.  We  believe,  then, 
that  there  exists  a set  of  state  variables  that  completely 
determine  profit  through  a fixed  real-world  relationship  W : 

v = W(d,j>) 

I call  such  a set  of  state  variables  "complete". 

We  believe  that,  if  we  could  predict  the  future  behavior  of 
the  state  variables  S_  and  if  we  knew  W , we  could  then  predict 
exactly.  We  are  thus  led  to  idealize  the  modeling  process. 

In  the  idealized  modeling  scheme,  we  determine  W by  examining 
objective  real-world  data  and  encode  our  uncertainty  about  the 
future  behavior  of  the  state  variables  as  a probability  distribution 
(sj^}.  Then,  through  the  relationship  W , this  uncertainty  about 
the  state  variable*  translates  directly  into  uncertainty  about 
profit,  expressed  as  the  profit  lottery: 


{v|d,4} 


<5(v  - W(d,S»  {S>} 


The  abridged  modeling  strategy  mimics  this  idealized  modeling 
scheme  in  that  no  uncertainty  is  attributed  to  the  dependence  relation- 
ship between  profit  and  the  state  variables.  Rather,  all  uncertainty 
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I I about  profit  is  attributed  solely  to  uncertainty  about  the  state 

variables  s acting  through  the  deterministic  model  g : 

" l <5(v  - g(d,s))  p(s) 

B_ 

But  this  arrangement  is  a misrepresentation  of  our  beliefs. 

The  modeling  that  we  actually  perform  in  a decision  analysis  falls 
short  of  the  Idealized  modeling  scheme  in  two  major  respects.  First, 
the  set  of  state  variables  s that  we  specify  for  use  in  the  model- 
ing strategy  is  generally  not  complete.  Typically,  in  a complex 
decision  problem,  although  we  can  identify  many  different  factors 
cuhL  influence  profit,  we  select  as  state  variables  only  those  few 
that  have  the  greatest  effect  on  profit  in  order  to  restrict  the 
assessment  and  computational  costs.  Consequently,  we  believe  that 
profit  depends  not  only  on  these  specified  state  variables  but  on 
other,  unspecified,  variables  as  well.  Even  if  we  knew  the  value  of 
each  specified  variable,  we  would  remain  uncertain  to  some  extent 
about  profit  becauue  of  the  effects  of  the  unspecified  variables. 

I 

Let  us  represent  the  unspecified  variables  as  the  vector  z. 

Then  si  and  z together  constitute  a complete  set  of  state  variables: 

(s,z)  ■ £ 

Assume  for  the  moment  that  we  know  W , the  real-world  re" ationship 
between  profit  and  this  complete  set  of  state  variables,  v ■ W(d,s,z). 
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Then,  our  uncertainty  about  the  dependence  relationship  between 
profit  v and  the  specified  state  variables  £ is: 

(v|d,s,4}  - [ {v|d,s,z,<S}  {z|s,4} 

L 

" | (v  - W(d,£,z>)  {z|s,4> 

Jz_ 

However,  we  generally  do  not  possess  sufficient  real-world 
u&ta  to  determine  W ; this  is  the  second  respect  in  which  we 
fall  short  of  the  idealized  modeling  scheme.  Although  we  believe 
that  there  does  exist  a fixed  real-world  relationship  between 
profit  and  the  complete  set  of  state  variables,  we  are  not  know- 
ledgeable enough  to  say  what  it  is. 

Let  W±  be  one  of  several  alternative  relationships.  Then, 
letting  (W± 1 4 > be  the  probability  that  W±  is  the  actual 

relationship,  our  uncertainty  due  to  this  insufficiency  of  data  can 
be  represented  as: 


{v|d,s,z,«S} 


J $(v  - Wi(d,s,£)){Wi|4} 


We  see,  then,  that  our  uncertainty  about  the  dependence  rela- 
tionship between  profit  v and  the  specified  state  variables  s 
has  two  sources  — our  uncertainty  about  the  effects  of  the 
unspecified  variables  z^  and  our  uncertainty  about  due  to  a 
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lack  of  sufficient  data: 


{v|d,8,4}  - [ [ 6(v  - W (d,s,z»  {w  |4>  {z|s,4} 

J1  1 

Clearly,  we  Ignore  both  of  these  sources  of  uncertainty  when 
we  assert  in  the  abridged  modeling  strategy  that  profit  is  com- 
pletely determined  by  the  specified  state  variables  jb  acting 
through  the  model  g : 

v - g(d,s) 

That  is,  in  making  this  assertion,  we  assume  clairvoyance  both  on 
£ given  £ and  on  W1  : 

{£|£,4}  ■ 6(z  - JlgCs))  for  some  fixed  Zq(b) 

and  ^ 1 4 * * 5 (w±  - Wk)  for  some  fixed  Wfc 

Assuming  this  clairvoyance,  we  would  have: 

{v|d,s,4}  - [ 6(v-  W (d,s,z»{W±|A}  {z  |s,4} 

' z_  1 

- | | 6(v  - Wi(d,s,z»  6(W1  - Wfc)  6(z  - z^s)) 

^ z 1 

“ - Wk(d,£,Zg(s))) 

- C(v-g(d,£)) 

where  g(d,s)  - Wfc(d,  8,^(3)) 
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But  we  do  not  possess  such  clairvoyance;  we  are  indeed 
uncertain  about  the  dependence  relationship  between  profit  and 
the  specified  state  variables.  Therefore,  we  should  regard  the 
approximation  made  in  the  abridged  modeling  strategy: 

{v|d,s,*}  = 6(v  - g(d,s)) 

as  just  that  — an  approximation  made  for  budgetary  reasons. 
Despite  the  similarity  in  appearances,  we  must  not  confuse  the  role 
of  the  model  g in  the  abridged  modeling  strategy  with  that  of  the 
relationship  W in  the  idealized  modeling  scheme.  Unlike  W , the 
model  g does  not  represent  a real-world  relationship  between 
profit  and  the  state  variables,  but  is  only  a convenient  surrogate 
for  the  probabilistic  relationship  thet  fully  expresses  our  uncer- 
tainty. The  model  g is  not  an  objective  description  but  rather  a 
subjective  expression. 

Several  writers  [l4, 17,26 J have  suggested  a probability- 
space  of  models  as  a vehicle  for  expressing  our  uncertainty  about  the 
dependence  relationship.  The  use  of  a model  space  is  attractive 
because  it  allows  us  to  employ  Bayesian  techniques  to  choose  among 
several  alternative  model  forms  in  the  face  of  real-world  data. 

If  we  use  such  a vehicle,  we  must  be  careful  how  we  interpret 
the  probability  assigned  to  each  model  in  the  space.  As  we  have  seen, 
unless  the  specified  set  of  state  variables  is  complete,  there  car  be 
no  "correct"  model  g that  deterministically  relates  profit  to  the 

state  variables;  therefore,  it  is  meaningless  to  speak  of  the 
probability  that  a particular  model  is  the  "correct"  one  or  that 
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it  embodies  the  actual  relationship. 

However,  there  is  a way  that  we  can  make  a meaningful  inter- 
pretation of  the  probability-space  of  models.  Suppose  that  we  are 
uncertain  about  which  of  several  relationships  w±  is  the  actual 
relationship  between  profit  and  the  complete  set  of  state  variables 
(s,z)  : 

{v|d,s,z,4}  =■  j 6(v  - Wi(d,£,£){Wi|4} 

Then,  we  can  write  the  conditional  mean  of  v given  d and  s as: 

* 

<v|d,s,4)  - v{v|d,s,4} 

v 

- J V f I 6(V  - W1(d,s,z)HWiUHz|s,4} 

V z 1 

■ [ f [ V 6(v  - W.  (d,s,z)){z|s,4} 

Ji  1 Jz  Jv  1 ““ 

= f tyj*}  j W (d,£,z){z.|£,A} 

= /±<W±U>  g/d.s) 

W.(d,£,z){zJjJ,4} 

L 

<v|d,s,WlP4> 


where  we  have  defined: 
gi(d,s) 


Thus,  we  can  expand  the  conditional  mean  <v|d,s,i)  over  a space  of 
models,  where  the  probability  assigned  to  model  gi  is  {W± ) ^ } , 

the  probability  that  the  actual  relationship  between  profit  and  the 
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complete  set  of  state  variables  is  , which  we  approximate 

by  gi- 

Note  that  even  if  W4  is  the  actual  relationship  between 
profit  and  the  complete  set  of  state  variables:  v = W^(d,s,z)  , 
it  is  not  true  that  gi  is  the  actual  relationship  between  profit 
and  the  subset  of  state  variables  £ : v * g^(d,s). 

To  say  that  there  is  no  "correct"  model  does  not  mean  that 
the  choice  of  function  g in  the  abridged  modeling  strategy  is 
unimportant.  The  choice  of  g partially  determines  the  closeness 
of  the  Dependence  Approximation: 


{v|d,s,4}  = 6{v  - g(d,s) ) 


and  thereby  affects  the  size  of  the  discrepancy  between  the  modeled 
and  authentic  profit  lotteries.  Although  no  function  makes  the 
approximation  exact,  we  want  one  that  makes  it  reasonably  close. 

An  obvious  candidate  is  the  conditional  expected  value: 


g(d,s)  = <v|d,s,<S> 

1.7  The  Adequacy  of  Modeling  in  Decision  Analysis 

We  perform  modeling  in  decision  analysis  not  for  its  own  sake, 
but  rather  to  help  us  with  the  analysis.  The  adequacy  of  the  model- 
ing in  a decision  analysis  must  therefore  be  judged  according  to  how 
well  it  serves  the  purposes  of  the  analysis. 
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In  a decision  analysis,  we  want  to  identify  the  optimal 
decision  — that  alternative  whose  corresponding  authentic 
profit  lottery  has  the  highest  expected  utility.  The  modeling  is 
adequate  if  it  leads  us  to  the  optimal  decision. 

It  is  easy  to  see  how  modeling  might  lead  us  astray.  The 
approximations  made  in  the  abridged  modeling  strategy  cause  us  to 
misstate  the  profit  lottery  and,  hence,  to  misstate  the  expected 
utility  corresponding  to  each  decision  alternative.  Hie  danger  thus 
exists  that  the  modeling  might  cause  us  to  mistakenly  identify  as 
optimal  an  alternative  whose  authentic  expected  utility  is  not 
highest. 

This  danger  is  especially  great  when  much  more  modeling  is 
required  to  obtain  the  profit  lottery  for  one  alternative  than  for 
another.  For  example,  suppose  that  the  two  decision  alternatives 
are  to  either  invest  or  not  invest  in  a particular  project.  Obtain- 
ing the  profit  lottery  for  the  first  alternative  may  require 
extensive  modeling  while  obtaining  the  profit  lottery  for  the  second 
may  require  none  (zero  profit).  In  this  case,  the  expected  utility 
for  the  first  alternative  will  be  misstated  while  the  expected 
utility  for  the  second  will  not,  and  our  decision  analysis  may  then 
identify  the  wrong  alternative  as  optimal.  (See  Figure  1.3). 

We  can  get  an  indication  of  the  adequacy  of  the  modeling  in  a 
decision  analysis  by  using  the  methodology  presented  in  Chapter  2. 
The  methodology  estimates  the  amount  by  which  we  misstate  the 


mm 


Figure  1.3s  Effect  of  Modeling  on  the  Decision 
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authentic  expected  utility  for  each  decision  alternative  because 
of  the  modeling  approximations.  If  the  misstatement  is  shown  to 
be  relatively  large,  we  may  presume  that  the  modeling  is  inade- 
quate . 

But  what  can  we  do  if  we  judge  the  modeling  to  be  inadequate? 
The  most  direct  approach  is  to  recalculate  the  profit  lottery  for 
each  alternative,  trying  to  get  it  '’closer”  to  the  authentic 
profit  lottery  by  tightening  one  or  both  of  the  modeling  approxi- 
mations: 


1.  {vjd,£,4} 

2.  {sJa} 


6(v  - g (d , s) ) (Dependence  Approx.) 

p(s)  (Distribution  Approx.) 


Hoping  to  minimize  the  additional  computational  costs,  we 
might  try  to  improve  on  these  approximations  without  enlarging  the 
decision  tree  by  choosing  a different  model  g in  the  first  approxi- 
mation or  by  choosing  different  discrete  values  of  the  state 
variables  in  the  second.  But  these  are  only  remedial  measures  and 
we  cannot  be  assured  that  they  will  get  us  closer  to  the  authentic 
profit  lottery.  Generally,  it  is  not  because  we  have  used  the 
"wrong"  model  or  the  "wrong"  discrete  values  of  the  state  variables 
that  the  modeling  is  Inadequate;  rather,  it  is  because  we  have  been 
forced  by  the  budgetary  constraints  into  making  tenuous  modeling 
approximations . 
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Hence,  to  get  closer  to  the  authentic  profit  lottery,  we  must 
enlarge  the  decision  tree.  We  can  make  the  tree  "taller"  by 
including  additional  state  variables  or  "bushier"  by  a 

finer  discretization  of  the  aleatory  variables  or  by  redesignat- 
ing some  of  the  non-aleatory  variables  as  aleatory.  Or,  acknowledg- 
ing that  we  are  uncertain  about  the  dependence  relationship,  we  can 
replace  the  deterministic  model  g with  a multi-point  discretiza- 
tion of  the  conditional  distribution  {v|d,s,4}.  Any  of  these 
measures,  of  course,  increases  the  computational  cost  of  obtaining 
the  profit  lottery. 

The  methodology  presented  in  Chapter  2 provides  us  with  an 
indirect  way  of  compensating  for  the  inadequacy  of  the  modeling. 

The  methodology  yields  an  estimated  correction  term  for  the  expected 
utility  of  each  alternative,  so  we  can  accept  as  optimal  that 
alternative  having  the  highest  corrected  expected  utility.  However, 
we  cannot  do  so  with  complete  confidence,  because  the  methodology  is 
not  exact. 

Clearly,  the  adequacy  of  the  modeling  in  a decision  analysis 
depends  not  only  on  our  modeling  skill,  but  also  on  the  computational 
resources  available  to  us.  Although  we  usually  like  to  think  that 
we  can  model  any  situation,  we  must  be  willing  to  concede  the  possi- 
bility that  budgetary  constraints  may  prevent  us  from  successfully 
modeling  a particularly  complex  decision  problem.  In  such  a case, 
the  results  of  the  decision  analysis  may  not  be  meaningful. 
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CHAPTER  2 


A METHODOLOGY:  THE  EFFECT  OF  MODELING  APPROXIMATIONS 

2 . 1 Introduction 

In  decision  analysis,  the  profit  lottery  that  we  compute  by 
modeling  is  not  identical  to  the  authentic  profit  lottery,  which 
we  ideally  want.  The  discrepancy  between  the  modeled  and  authentic 
profit  lotteries  is  due  to  approximations  made  in  the  modeling 
process.  The  methodology  developed  here  allows  us  to  quantita- 
tively relate  the  size  of  the  discrepancy  to  the  modeling  approxi- 
mations made  in  the  analysis. 

As  an  overriding  goal,  I have  tried  to  keep  the  methodology 
simple  enough  that  it  can  be  employed  quickly  and  easily  to  check 
the  adequacy  of  the  modeling  in  any  decision  analysis. 

Consider  a decision  problem  defined  by  decision  variable  d 
and  outcome  variable  v.  We  would  like  to  obtain  the  authentic 
profit  lottery  {v|d,*}.  To  do  this,  we  first  specify  a set  of 
state  variables  b on  which  we  believe  v partially  depends.  The 
dependence  of  v on  d and  s.  can  be  expressed  as  a conditional 

probability  distribution  {v|d,s_,4>. 

If  we  could  assess  the  authentic  distributions  {sj*}  and 
{v|d,s_,A>  and  if  we  could  perform  the  necessary  integration,  we 
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would  obtain  the  authentic  profit  lottery  via  the  expansion 

j 

equation: 

> 

(2.1)  {v|d,4}  = {v|d,s,4}{s|4} 

S 

— 

However,  this  is  not  generally  possible.  Usually,  we 
cannot  afford  the  tine  and  effort  to  fully  assess  the  distri- 
butions nor  can  we  analytically  perform  the  required  integration 
over  arbitrary  continuous  functions. 

Consequently,  we  make  two  modeling  approximations  to  simplify 
our  task.  First,  we  assume  that  the  probability  distribution  on  £ 
can  be  represented  by  a discrete  probability  function  p(s).  I call 
this  the  Distribution  Approximation: 

(2.2)  {sJa } = p(s) 

Included  in  the  Distribution  Approximation  is  the  fixation  of  each 
non-aleatory  variable  at  a single  representative  value,  in  effect 
approximating  its  distribution  by  a single-point  discretization. 

Secondly,  we  assume  that  the  dependence  of  v on  d and  _s 
can  be  represented  by  a deterministic  function  g.  I call  this  the 
Dependence  Approximation: 
v * g(d,s) 
jr,  equivalently, 

(2.3)  {v|d,£,4}  = 6(v  - g(d,s)) 
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Having  made  these  two  modeling  approximations,  we  obtain  the 
familiar  block  diagram  representation  of  the  decision  problem  and, 
equivalently,  the  decision  tree,  as  shown  in  Figure  2.1.  We  can 
then  easily  compute  the  modeled  profit  lottery: 

(2.4)  {v|d,4)m  - £ fi(v  - g(d,s) ) p(s) 

8 


But,  because  of  the  modeling  approximations  (2.2)  and  (2.3)  , 

the  modeled  profit  lottery  {v|d,4}  from  (2.4)  is  not  identical 
to  the  authentic  profit  lottery  {v|d,4}  from  (2.1). 

In  the  methodology  that  follows,  we  shall  quantitatively 
characterize  the  modeling  approximations  and  determine  their  effects 
on  the  discrepancy  between  the  profit  lotteries.  We  shall  look  at 
the  approximations  one  at  a time,  starting  with  the  Distribution 
Approximation. 

To  facilitate  the  discussion,  I shall  use  subscripts  for  the 
outcome  variable  v to  denote  the  various  stages  of  the  develop- 
ment of  the  methodology  (see  Figure  2.2): 


v *=  the  modeled  outcome  variable 
m 

v ■ the  outcome  variable  corrected  only  for 
x 

the  Distribution  Approximation 

v ■ the  outcome  variable  corrected  for  both 
a 

approximations  (i.e.,  the  authentic  outcome 
variable,  defined  as  the  variable  whose  distribution 
is  the  authentic  profit  lottery) 
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v - g(dlts3) 


Decision  Tree 


Modeled  outcome  variable 


v * 
m 

fv  |d,A>  “ I<S(v  - g(d,s) ) p(£) 
s 


* 

Correction  for  the 
Distribution  Approximation: 

{s  | A } - p (s) 


V 

x 


Outcome  variable  corrected  only 
for  the  Distribution  Approximation 


* 

6(v  - g(d,s)){sj.6} 
s 


i 

Correction  for  the 
Dependence  Approximation: 

{v|d,£,4>  a 6(v  - g(d,s)) 


♦ 


v ■ Authentic  outcome  variable 
a 

■Cv  | > = [ {v|d,s,4}{£|4> 


The  goal  of  the  methodology  is  to  determine  correction  terms 
Av  and  Av  for  the  moan  and  variance  of  the  outcome  variables 

Av  ■ v - v 

a m 

v V v 

Av  = v - v 

a m 

We  shall  consider  In  Sections  2.2,  2.3  and  2.4  the  case  In 
which  there  Is  only  one  state  variable  s and  then,  In  Section  2.5, 
we  shall  extend  the  results  to  the  multivariate  case. 

2.2  Effect  of  the  Distribution  Approximation 

Assume  that  the  state  variable  8 is  scalar  and  assume  for  the 
moment  that  v * g(d,s)  exactly.  We  want  to  determine  the  effect 
on  the  profit  lottery  of  using  the  discrete  distribution  p(s)  on 
the  state  variable  rather  than  the  authentic  distribution  {s | A } . 

For  notatiorial  ease,  let 

n = the  state  variable  whose  distribution  is  p(s) 
x ■ the  state  variable  whose  distribution  is  {s|4} 

The  choice  of  "n"  and  "x"  emphasizes  that  p(s)  is  ulscrete 
while  { s 1 4 } is  generally  continuous. 

We  characterize  the  Distribution  Approximation  by  the  differences 
in  the  mean  and  variance  of  x and  n : 

As  ■ x - n 

V V V 

As  m x - n 
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(2.9)  vx  = g(d,n)  + g'  (d,n) As  + V.g"(d,n)(x  + As2) 

(2.10)  vx  = g’(d,n)2  x + ^"(d.n)2  (x2  + 2xAs2) 

+ 2g * (d,n)g"(d,n)  xAs 


By  subtraction,  we  get: 


(2.11) 

✓-N 

< i 

X 

i 

0*' 

'w' 

= g'(d,n>As  + *5g"  (d,n)  (As  + As2) 

(2.12) 

/—V 

>>B 

1 

X 

> > 

- g'(d,n)2As  + 1sg"(d,n)2j^As(x+n)  + 2xAs 

+ 2g* (d,n)g"(d,n)  xAs 

These  equations  yield  the  estimated  correction  terms  for  the 
effects  of  the  Distribution  Approximation.  Equation  (2.11)  is 
exact  for  quadratic  g and  (2.12)  is  exact  for  linear  g. 

2.3  Effect  of  the  Dependence  Approximation 

We  now  consider  the  effect  on  the  profit  lottery  of  the 
Dependence  Approximation.  To  compute  the  profit  lottery  {v|d,4}m  , 

we  assume  that  v is  completely  determined  by  d anu  s : 

v “ g(d,s) 

In  most  real  situations,  however,  we  believe  that  v also  depends 
on  other  less  important,  vinspecified  variables.  Hence,  by  assuming 
the  deterministic  relationship,  we  have  suppressed  whatever  uncer- 
tainty wt;  may  have  about  v for  fixed  values  of  d and  s. 


40 


We  characterize  the  Dependence  Approximation,  then,  by  a 
random  variable  e : 


v - g(d,s)  + e 

Our  residual  uncertainty  about  v given  d and  s can  be 
expressed  as  a conditional  probability  distribution  on  e , 

'.eid, s, 4}.  For  the  purposes  of  the  methodology,  we  assess  the 
conditional  mean  (e|d,s,4)  and  variance  \ e|d,s,4). 

Recalling  that  vx  is  corrected  for  the  effects  of  the  Dis- 
tribution Approximation,  we  have: 

(2.13)  va  - g(d,x)  + e - vx  + e 

Taking  the  mean  and  variance  of  this  equation  (see  Appendix  A.  3), 
we  get: 

(2.14)  (v  - v ) - e 

u X 

(2.15)  (v  - v ) “ e + 2cov(v  ,e) 

ax  x 

These  equations  yield  the  correction  terms  for  the  effects  of 


the  Dependence  Approximation. 
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We  can  calculate  e , e and  cov(v  >e)  as  follows  (see 


Appendix  A. 4)  : 

(2.16)  l - ja<e|d,s,4>{s|4} 


(2.17) 


v 

e « 


-2 

e 


|v<e|d,s,4)  + <e|d,s,4)2  J {s|4>  - 

^s  L -* 

cov(vx>e)  ■ |s  |g(d,s)  <e|d,s,i>  {s|4> j - vx< 


(2.18) 

Note  that  each  of  these  calculations  requires  the  integration 
of  a function  of  d and  s (e.g.,  (e|d,s,4))  over  the  authentic 

distribution  {s | •*} . In  most  cases,  we  are  unable  to  perform  this 
integration  and  must  approximate  it  with  a summation  over  the  dis- 
crete distribution  p(s).  This  is  just  another  utilization  of  the 
Distribution  Approximation;  that  is,  for  an  arbitrary  function 
f(d,s)  , 

t 

f(d,8)  (s|a>  - \t (d,s)  p(s) 


We  can  use  the  results  of  Section  2.2  to  find  correction  terms  for 
these  calculations.  Letting: 


x 


f (d,x)  and  53  f(d,n) 


we  see  that 


x 


f(d,s)  {s|4> 


l f(d,s)  p (s) 
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Then,  we  expand  z and  z about  n , take  the  mean  and  sub- 
’ r x n 

tract: 

(2.19)  (i  - zn)  a f’(d,n)As  + 4f"(d,n)(As  + As2) 

This  equation,  with  appropriate  substitution  for  the  function 
f (d,s)  , yields  the  estimated  correction  term  for  each  of  the 

calculations  in  (2.16)  , (2.17)  and  (2.18)  to  compensate  for 

the  use  of  p(s)  rather  than  {s|a}. 

Two  special  cases  of  the  Dependence  Approximation  merit  our 
attention.  As  the  first  special  case,  suppose  that  e is  indepen- 
dent of  s ; that  is,  suppose  that  our  residual  uncertainty  about 

v given  d and  s does  not  depend  on  s.  Then,  we  need  only 
- v 

assess  e and  e to  characterize  the  Dependence  Approximation: 


<e|d,s,4>  - 

<e|d,4 ) - <e 

'<e|d,s,-&>  - 

v , V 

<e|d,4>  = e 

If  e is  independent  of  s , it  is  also  independent  of  vx  , so: 
cov(vx,e)  - 0 

Thus,  when  e is  independent  of  s , the  correction  terms  for  the 
Dependence  Approximation  from  (2.14)  and  (2.15)  become: 
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(2.20)  <vfl  - vx)  - e 

V V V 

(2.21)  (va  - vx)  - e 


As  the  second  special  case,  suppose  that  we  find  it  easier  to 
assess  our  residual  uncertainty  about  the  model  g in  terms  of  its 
output  vx  rather  than  its  input  s.  Then,  we  would  assess  <e|d,vx,i> 
and  V<(e|d,vx,4)  instead  of  (e|d,s,A^  and  (e|d,s,4)  . In  this 
case,  we  can  always  determine  (e|d,s,<i)  and  (e|d,s,A)  by  a 
simple  change  of  argument: 

<e|d,s,4>  - <e|d,vx,4> 

V<e|d,s,4>  - V<e|d>vx«4)  v -g(d,s) 

x 

However,  suppose  that  our  residual  uncertainty  about  the  model 
is  proportional  to  its  output  vx  in  the  following  sense . 

<e|d,vx,4>  = A vx 

V<e|d.vx.^>  " B vx2 

where  A and  B are  constants.  Such  would  be  the  case,  for  example, 
if  {e|d,vx,4}  were  uniform  or  triangular  over  e with  a base  pro- 

portional to  v . Then,  we  have  (see  Appendix  A. 5)  : 


<e|d,4> 


v v 
e “ > 


'<e  |d,<S> 


2 v - 2 

(B  + A ) vx  + B vx 


cov(vx,e) 


So,  for  this  special  case,  the  correction  terms  for  the  effects  of 
the  Dependence  Approximation  from  (2.14)  and  (2.15)  become: 


(2.22)  (va  - vx) 


A v 

x 


(2.23)  (vfl  - vx) 


(B  + 2A  + A2) 


v 

V 

X 


+ 


2.4  Combined  Effect  of  Both  Modeling  Approximations 

- v 

The  correction  terms  Av  and  Av  for  the  combined  effect  of 
both  modeling  approximations  are  simply  the  sums  of  the  correction 
terms  for  each  approximation: 


(2.24)  Av 


(v  - v ) 
a m 


(v  - v ) 
a x 


(v  - V ) 
x m 


(2.25)  Av  = 


v V 

(v  - V ) 
a m 


V V 

(v  - V ) 
a x7 


Effect  of 
Dependence 
Approximation 
(2.14),  (2.15) 


V V 

(v  - V ) 
x nr 


Effect  of 
Distribution 
Approximation 
(2.11),  (2.12) 


In  most  cases,  the  deterministic  model  that  we  use  to  compute  the 
profit  lottery  is  composed  of  several  submodels.  We  can  use  the  methodol- 
ogy sequentially  for  each  submodel  to  determine  the  cumulative  effect  on 
the  profit  lottery  of  the  modeling  approximations  made  in  all  of  the  sub- 
models. For  example,  referring  to  the  model  shown  in  Figure  2.3,  we 

would  first  use  the  methodology  to  determine  the  correction  terms  Ay 
v 

and  Ay  for  the  intermediate  variable  y.  Then,  considering  y as  an 
input  variable  to  the  second  submodel,  we  would  use  the  methodology  again 
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to  determine  the  correction  terms  Av  and  Av  for  the  outcome 
variable  v. 

Note  that,  if  we  use  the  methodology  sequentially  for  several  sub- 
models, we  need  the  modeled  mean  and  variance  for  each  intermediate 
variable.  Generally,  the  computer  programs  used  to  calculate  the 
modeled  profit  lottery  can  easily  be  arranged  to  report  these  quantities. 

Going  beyond  the  profit  lottery,  we  can  determine  the  effect  of  the 
modeling  approximations  on  the  expected  utility  and  the  certain  equiva- 
lent by  considering  the  utility  function  u(v)  as  just  another  submodel. 
We  want  to  find  the  effect  on  the  expected  utility  u of  using  the 
modeled  profit  lottery  (v|d,4}n  rather  than  the  authentic  profit 
lottery  {v|d,4}  as  the  distribution  on  input  variable  v.  This  is 
the  Distribution  Approximation  once  again,  which  we  characterize  by  the 

__  V 

correction  terms  Av  and  Av  . From  (2.11)  , we  have: 

(2.26)  Au  - u' (v  )Av  + W (v  ) (Av  + Av^) 

m ra 

Suppose,  for  example,  that  the  utility  function  is  exponential 
with  risk  coefficient  Y : 

- Yv 

u(v)  = 1 - e 


Then: 


u'  (v) 
u"(v) 
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Substituting  into  (2.26)  , we  have: 

Au  = ye  ^Vm(Av)  - e (Av  +Av  ) 

= Ye”^Vm^Av  - **y(Av  + Av^)^j 

The  corrected  expected  utility  is  then: 

u * u + Au 
a m 

The  corrected  certain  equivalent  can  be  computed  from  the  corrected 
expected  utility. 

2.5  Extension  to  Multivariate  Models 

We  can  extend  the  methodology  developed  in  the  previous  sections 
to  decision  problems  in  which  there  are  more  than  one  state  variable. 
However,  the  methodology  is  practical  only  when  the  state  variables  are 

mutually  independent. 

Consider  the  multivariate  model  of  a decision  problem  shovn  in 
Figure  2.4.  Outcome  variable  v is  dependent  on  decision  variable  d 
and  on  the  state  variable  vector  s.  The  dependence  is  fully  expressed 
by  the  conditional  probability  distribution  {v|d,s,4)m,  but  is 
approximated  for  purposes  of  analysis  by  a deterministic  relationship 

g * 

{v|d,£,A>  - 6(v  - g(d,s.)) 

Also,  the  authentic  distribution  on  s.  , {s^l-i}  , is  approximated  by 

a discrete  probability  function  p(s). 


Having  made  these  two  modeling  approximations,  we  compute  the 
modeled  profit  lottery: 

{v|d,4}m  - Jfi(v  - g(d,s))  p (s_) 


We  want  to  determine  the  discrepancy  between  this  modeled  profit 
lottery  and  the  authentic  profit  lottery: 

{v|d,4}  * J {v|d,£,.&}  {sj.6} 

f a 


We  develop  the  methodology  as  in  the  previous  sections,  looking 
first  at  the  Distribution  Approximation  and  then  at  the  Dependence 
Approximation. 


n - the  state  variable  vector  whose  distribution 
is  p(s) 


x ■ the  state  variable  vector  whose  distribution 


is  Q.1  -A  ) 


We  characterize  the  Distribution  Approximation  by  the  differences 
in  the  mean  and  variance  of  x and  n : 

As  ■ x “ E 

v V V 

As  ■ x - n 


Then,  we  expand  both  vffl  * g(d,n)  and  vx  ■ g(d,x)  about  n and 
take  the  mean  and  variance  (see  Appendix  A. 6),  using  simplifying 
assumptions  as  before  to  eliminate  third-  and  fourth-order  terms  (see 


50 


Appendix  A. 2).  As  shown  in  Appendix  A. 6,  if  the  state  variables  are 
not  mutually  independent,  the  resulting  expressions  are  not  useful 
because  they  require  more  information  than  we  have  about  the  distri- 
butions on  s_.  So,  assuming  that  the  state  variables  are  mutually 
independent , we  have : 


<vx  “ vm)  ~ Igi(d,n)Asi  + ^ ^g11(d,n)Asi 


(2.27) 


(2.28) 


where 


+ %l  1*^4  (d,n) As^Ae. 

ij 


V V r - 2 v 

(vx  " Vm)  = ASi<d*H>  Asi 
i 


_ _2VV  v v v _ 2 v _ 2 

+ h\  ).g1;j(d,n)  (x^-n^j+x^Sj  +XjAsi  ) 

+ III  gtj (d,n)glk(d,n)  x±As  Ask 

i i k*j 

- - V — 

+ 2 4 Zg±  <d »IL) g± j (d ,n>  x1As. 

1 J 


g1(d,n)  - — g(d,.s) 


s ■ n 


Eij<d’i>  ' s,V8..  8<d>-) 


i J 


s = n 


Equations  (2.27)  and  (2.28)  yield  the  estimated  correction 


terms  for  the  effects  of  the  Distribution  Approximation. 


As  e.  op  zial  case,  suppose  that  there  are  just  two  Independent 


state  variables  s.  and  e..  Then  th-.  correction  terms  from 
(2.27)  and  (2.28)  reduce  to: 


(vx  " \)  ~ g1(d,n)As1 

+ *5g13  (d,n)(As1+As12) 


s 


1 terms 


(2.29) 


+ g2(d,n)As2 


+ Jsg22(d,n)(As2  +As22) 


s 


2 terms 


+ JSg12(d,n)As1^2 


crosb  term 


v v 

(v  - v ) 
v x m' 


- 2 v 

g1(d,n)  Asx 


_2lvvv  v _ 2~| 

•*  *5gi:L(d»n)  ^^(x^j^)  + 2x1As1  j 
+ 2g1(d,n)g11(d,n)  x^^ 


S1 

terms 


- 2 v 
+ g2(d,n)  As2 

+ Hh2 


Cd,n)2|is 


V V 
2(x2+n2)  + 2x2 


+ 2g2(d,n)g22(d,n)  x2As2  j 


(Expression  continued  on  next  page) 
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(2.30) 


— 2 v v vv  v _ o v 2 
+ g12(d,n)  (x^  - n1n2‘hc1^2  +*2^1  5 

+ 2811(d»i)g12(d»!L)  x1As1As2 

+ 2g22(d.n)g12(d,n)  x^s^A^ 

+ 2g1(d,n)g12(d,n)  x1As2 

- _ v _ 

+ 2g2 (d,n)g12 (d,n)  x2As1 

Since  the  outcome  variable  is  still  a scalar,  we  handle  the 
Dependence  Approximation  exactly  as  in  Section  2.3.  We  characterize 
the  approximation  by  a scalar  random  variable  e : 

v - g(d,s)  + e 

and  assess  the  conditional  mean  and  variance  (e[d,s,4)  and 
v<(e|d,£,4).  Then,  the  correction  terms  for  the  effects  of  the 
Dependence  Approximation  from  (2.14)  and  (2.15)  are: 

(2.31)  (v  - v ) - e 

a x 

(2.32)  (v^  - v^)  ■ e + 2cov(vx>e) 

_ v 

where  e , e and  cov(vx>e)  are  calculated  as  in  (2.16)  , 

(2.17)  and  (2.18). 

If  we  use  the  discrete  distribution  p(s)  instead  of  the 
authentic  distribution  { s_| A )■  in  the  calculations  of  e , e and 
C0v(vx*e>  * We  can  determine  correction  terms  for  these  calculations 
as  in  (2.19).  The  multivariate  counterpart  of  (2.19)  is: 


cross 

terms 


IUIUWPIIWIU^IIIU.1 


'.W’<  " 


I 1. 11,4.1-1  LU 


C\-zn)  = l f <d,n)AS  + h l f (d,n)As 
i i 1 

(2* 33)  + H l l f..(d,n)As  As. 

i j 2 13 

Finally,  the  correction  terms  for  the  combined  effects  of 


both  modeling  approximations 

are: 

Av  - (v  - v ) = 

a m 

(v  - v ) + 

a x 

/'-N 

■>s 

1 

1 >* 
Vw/ 

V V V 

Av  * (v  - v ) - 

a m 

V V 

(v  - V ) + 

a x 

V V 

(v  - V ) 
x m 

Effect  of 

Effect  of 

Dependence 

Distribution 

Approximation 

Approximation 

(2. 31), (2. 32) 

(2. 27), (2. 28) 
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AN  EXAMPLE:  COMPETITIVE  PRICING  DECISION 

3.1  Introduction 

The  decision  making  client  plans  to  produce  and  sell  a new  home 
appliance  during  the  coming  year.  He  knows  that  he  will  have  one  major 
competitor  and  that,  because  the  two  competing  appliances  are  virtually 
Identical  In  function,  the  share  of  the  market  each  competitor  captures 
will  be  determined  primarily  by  the  relative  selling  prices.  The 
client  is  not  sure  what  price  his  competitor  will  set. 

Furthermore,  he  is  uncertain  about  the  size  of  the  total  market. 

He  believes,  however,  that  it  will  depend  strongly  on  the  lower  of  the 
two  selling  prices  and  on  the  general  economic  conditions  during  the 
year. 

Finally,  uncertainty  also  surrounds  the  manufacturing  costs.  The 
appliance  contains  an  expensive  raw  material,  but  the  client  is  uncertain 
how  much  the  material  will  cost  him  and  how  much  of  the  material  is 
required  for  each  unit. 

The  client's  decision  is  which  selling  price  he  should  set  for  his 
appliance.  Section.  3.2  is  a description  of  the  decision  analysis  of 
this  problem  and  Section  3.3  is  an  application  of  the  methodology  of 
Chapter  2 to  determine  how  much  effect  the  modeling  approximations  made 
in  the  analysis  have  on  its  results. 
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3.2  Decision  Analysis  of  the  Problem 


Deterministic  phase:  In  the  analysis  of  this  problem,  we  first 

construct  a deterministic  model,  as  shown  in  Figure  3.1.  The 
decision  variable  is  the  selling  price  p and  the  outcome  variable 
is  the  profit  s.  Given  p and  the  competitor's  selling  price  y , 
the  client's  share  of  the  market  z is  determined  by  the  Market  Split 
Model : 


P 1 y 

"h“e  k ■ mTsT  * 5'68 

p < y 


See  Figure  3.2  for  a graph  of  the  Market  Split  Model. 

The  Demand  Model  determines  the  total  market  size  n given  the 
lower  selling  price  x and  an  indicator  of  the  general  economic 
conditions  f : 


n 


(millions  of  units) 


See  Figure  3.3  for  a graph  of  the  Demand  Model. 
The  client 's  total  sales  q is  then: 


3 “ z n (millions  of  units) 

His  total  revenue  r is : 


r - p q 


(millions  of  dollars) 
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Figure  3.2:  Market  Split  Iladel 


z = market  share 
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The  manufacturing  cost  per  unit  v depends  on  the  cost  h of 
the  primary  raw  material,  on  the  per  unit  usage  factor  w of  that 
material,  on  other  variable  costs  b and,  because  of  a volume  dis- 
count on  the  primary  raw  material,  on  total  sales  q : 

v = b + (.85  + .15eKqh)  h w (dollars) 

where  K = -In (81)  - -4.4 

Then,  given  fixed  costs  t , the  total  cost  c is: 

C - t + V q (millions  of  dollars) 

Finally,  profit  s is  simply: 

s = r “ c (millions  of  dollars) 

Next,  we  assign  low,  nominal  and  high  values  for  the  state  vari- 
ables, as  follows: 


State  variable 

Low 

Nominal 

200 

U-f 

Competitor's  price: 

y ($) 

180 

nlgn 

250 

Economic  indicator: 

f 

0.75 

1.00 

1.25 

Material  usage: 

w 

1.10 

1.25 

1.40 

Raw  materia]  cost: 

h ($) 

30 

40 

50 

Other  var.  costs: 

b ($) 

30 

40 

50 

Fixed  costs: 

t ($106) 

7 

9 

11 

I 


f 

j 


5 

l 

i 


j 
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Using  the  deterministic  model  and  the  nominal  values  of  the  state 
variables,  we  compute  the  deterministic  profit  for  each  of  several 
different  values  of  the  decision  variable  p.  The  results  are  shown 
on  Figure  3.4  The  highest  deterministic  profit  of  s = $6.97  million 
corresponds  to  the  decision  p ■ 170. 

We  next  calculate  the  deterministic  sensitivity  of  the  profit  to 
each  of  the  state  variables  by  varying  each  state  variable  through  its 
range  while  holding  the  others  at  their  nominal  values.  The  results  of 
the  deterministic  sensitivity  are  shown  in  Table  3.1.  We  find  that  pro- 
fit is  relatively  insensitive  to  the  material  usage  factor  w and  to 
the  fixed  costs  t.  Profit  is  highly  sensitive  to  the  economic  indicator 
f , although  the  deterministically  optimal  decision  is  not  affected  by 
f , as  it  is  by  the  competitor’s  price  y , the  raw  material  cost  h 
and  the  other  variable  costs  b. 

Probabilistic  phase;  Because  deterministic  profit  is  relatively 
insensitive  to  changes  in  w and  t , we  designate  thorn  as  non-aleatory 
and  fix  them  at  their  nominal  values.  We  designate  ’"he  remaining  state 
variables  as  aleatory  variables  and,  in  interviews  with  the  client, 
assess  the  probability  distribution  for  each  aleatory  variable,  as  shown 
in  Figures  3.5  through  3.8.  We  then  approximate  each  of  these  distri- 
butions by  a discrete  distribution  as  shown  in  the  figures  to  allow 

assignment  of  probabilities  to  the  decision  tree.  The  competitor’s  price 

y is  discretized  into  eight  steps,  the  economic  indicator  f into  five 

steps  and  the  raw  material  cost  h and  the  other  variable  costs  b into 

three  steps  each.  The  resulting  decision  tree  has  360  "tips”  for  each 
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Authentic  distr 


Discrete  approx 


Figure  3.? 


Probability  Distribution  on  b 


b - other  variable  costs 


decision  alternative. 


Computations  from  the  decision  tree  result  in  the  modeled  profit 
lottery  {s|p,4}^  for  each  of  several  different  values  of  the 
decision  variable  p.  We  determine  the  relative  desirability  of  each 
profit  lottery  by  means  of  the  utility  function.  Interviews  with  the 
client  reveal  that  his  utility  function  on  profit  s is  exponential 
with  risk  aversion  coefficient  .1  : 

u(s)  = 1 - e 8 for  s in  millions  of  dollars 

Using  the  utility  function,  we  compute  the  expected  utility  and  the 
certain  equivalent  of  each  profit  lottery : 

- l {s|p,4}mu(s) 

8 

CEm(p)  = u_1(up)  - -10  ln(l  - Gp) 

Figure  3.9  is  a graph  of  the  certain  equivalents.  The  most  preferred 
profit  lottery,  shown  in  Figure  3.10,  corresponds  to  p = 160  and  has 
a certain  equivalent  of  CE  - $6.47  million. 


Informational  phase:  We  calculate  the  probabilistic  sensitivity 

of  the  certain  equivalent  to  each  aleatory  variable  by  holding  the 
variable  to  each  of  its  values  while  allowing  the  other  aleatory  variables 
to  remain  probabilistic.  Table  3.2  chows  the  results  of  this  probabi- 
listic sensitivity  analysis.  As  in  the  deterministic  sensitivity 
analysis,  the  optimal  decision  is  not  affected  by  th^  economic  indicator 
f.  The  optimal  decision  is  most  sensitive  to  the  competitor's  price  y. 
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Using  the  results  of  the  probabilistic  sensitivity  analysis,  we 
compute  the  value  of  clairvoyance  for  ecch  aleatory  variaole,  as 
follows: 


Aleatory  variable 
Competitor’s  price:  y 

Economic  indicator:  f 

Raw  material  cost:  h 

Other  var.  costs:  b 


Value  of  clairvoyance 
$378,000 
0 

$ 46,000 
$ 31,000 


3.3  Effect  of  the  Modeling  Approximations 

We  can  now  consider  what  effect  the  modeling  approximations  made 
in  the  analysis  have  on  its  results.  The  fundamental  question  we  ask 
here  is  how  wel1  each  of  the  profit  lotteries  computed  in  the  analysis 
(e„g..  Figure  3.10)  expresses  the  client's  beliefs  about  the  future 
consequences  of  his  decision  or,  equivalently,  how  close  the  computed 
profit  lottery  is  to  the  authentic  profit  lottery. 

The  discrepancy  betveen  the  computed  and  the  autnentic  profit 
lotteries  arises  from  two  types  of  modeling  approximations:  1)  the 

approximation  of  the  probability  distribution  on  each  state  variable  by 
a discrete  distribution  (Distribution  Approximation)  and  2)  the 
approximate  of  the  probabilistic  relationship  between  variables  by  a 
deterministic  relationship  (Dependence  Approximation). 

We  can  use  the  methodology  developed  in  Chapter  2 to  estimate  the 
discrepancy  in  the  mean  and  variance  of  the  outcome  variable,  profit, 
caused  by  these  approximations. 
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We  characterize  the  Distribution  Approximation  for  each  state 
variable  by  the  differences  in  the  mean  sr.d  variance  of  the  variable 
using  the  two  distributions,  as  shown  in  Figures  3.5  through  3.8. 

These  correction  terms  for  all  of  the  state  variables  are  shown  in 
Table  3.3. 

We  shall  trace  the  effects  of  the  various  modeling  approximations 
through  the  deterministic  model  shown  in  Figure  3.1,  using  the  metho- 
dology sequentially  for  each  submodel.  Figure  3.11  shows  the  flow  of 
this  sequential  use. 


Effect  on  market  share  z of  the  Distribution  Approximation  on  competi 
tor’s  price  y : 

From  the  Market  Split  Model,  we  have  for  p £ y : 

z - 1 - where  k - 

Taking  derivatives  with  respect  to  y : 


z’(y)  - 


*”<y>  - 


-(k+2) 


For  p » 160  , these  derivatives,  evaluated  at  y ■ 201,1  are; 

m 


z’  (y  ) = 38.54  x 10' 

wm 

z"(y  ) = -1.28  x 10' 


•feble  3.3.  Correction  Terms  for  Distribution  Approximations 


Authentic 

Modeled 

Difference 

Competitor1  s 

y - 201.10 
a 

y = 201.10 
m 

Ay  = 0 

price 

V 

V 

V 

ya  - 310.46 

d 

y = 293.80 
m 

Ay  = 16.66 

Economic 

t = 1.003 

cL 

f = 1.017 
m 

Af  = -0.014 

indicator 

V 

V 

V 

f = 0.025 
a 

f = 0.020 

m 

Af  = 0.005 

Raw  material 

h - 39.29 
a 

h = 39.10 
m 

Ah  = 0.19 

cost 

V 

V 

V 

h » 57.99 
a 

h = 54.19 
m 

Ah  = 3.80 

Other 

Sa  - 41.33 

b = 40.48 

Til 

Ab  = 0.85 

variable 

V 

V 

V 

costs 

b = 32.89 
a 

b = 31.77 
m 

Ab  = 1.12 

Material 

wa  - 1,25 
a 

w = 1.25 
m 

Aw  = 0 

usage 

V 

V 

V 

factor 

w = 0.01 
a 

w = 0 
m 

Aw  = 0.01 

Fixed 

t - 9.0 
a 

t = 9.0 
m 

At  = 

costs 

V 

V 

V 

t = 1.04 
a 

t = 0 
m 

At  = 1.04 
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Cumulative 
effect  on 
Certain 

equivalent  CE 
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Then,  noting  that  Ay  = 0 , we  have  from  Equations  (2.11)  and 

(2.12)  the  correction  terms  for  the  effect  on  z of  the  Distribu- 
tion Approximation: 


<zx  - *„>  - h 

(*X  ' ZJ  * z'(>'»)2  £ * H 2"v2  4y(ya4m) 

V V V 

Substituting  from  Table  3.3  the  values  of  y&  , y and  Ay 


" im)  = ( 8*33)  z"(ym) 


(zx  ~ V = <16*66>  z’(yJ2  + (5033.5)  z"(y  ) 


For  p = 160  , these  correction  terms  are: 


<**-*„>  - --“Oil 


V V _A 

(z  - z ) = 3.30  x 10  H 

x m 


Effect  on  market  share  z of  the  Dependence  Approximation  in  the 


Market  Split  Model: 


The  model  asserts  that  the  spilt  of  the  market  is  determined 
solely  by  the  two  selling  prices.  The  client  believes,  however,  that 
the  model  does  not  fully  represent  his  uncertainty  about  the  market 
^or  instance,  he  believes  that  other  factors,  such  as  adver- 
tising and  brand  loycaty  may  also  influence  the  market  split,  although 
to  a lesser  extent  than  the  relative  price?. 
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We  characterize  this  Dependence  Approximation  by  the  random 


variable  e : 
z 


z = z + e 
a x z 


The  client  believes  that  the  model  is  unbiased  but  that  his  residual 
uncertainty  has  a standard  deviation  of  about  8%  of  the  value  of  z 
given  by  the  model  for  z <_%  and  a symmetric  amount  for  z > %. 
Therefore,  he  assesses: 


<e_  | P » z ,4)  = 0 

z x 


V<ez|P.zx.4> 


1 2 
z 


150  x 


0 < z < h 
— x — 


V 


1 ,2  , 

150  (1  " 2X*  * - zx  - 1 


Since  the  conditional  mean  is  zero,  e » 0 and  cov(z  ,e  ) = 0 

z x z 


v 

Then,  calculating  ez  for  z > % : 


v 

e 


/<ezlp«zx’4>  {zx!p«4> 


’ I Bo  a - V2f\lp.i) 


150  (1  - 2z  + z 2) 

X X ' 


TT?T  (1  - 2z  + z + z 2) 
150  x x x ' 


From  (2.14)  and  (2.15)  , the  correction  terms  for  the  effect  on 

z of  the  Dependence  Approximation  are: 


(‘a  " ^ 


za  “ ZJ  “ e„  + 2cov(z  ,e  ) 


X z 


V 

e 
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For  p ■ 160 


the  modeled  mean  and  variance  of  z are : 


z - . 847 

m 

z =*  38.45  x 10" 

m 

So,  for  p = 160  : 


z - z + (z  - z)  = .846 

x in  x in 


V V V V _/. 

z = z + (z  - z ) ■ 41.75  x 10 


x m x m 


Then,  the  correction  terms  for  the  efilect  on  z of  the  Dependence 
Approximation  are  for  p ■ 160  : 


(z  z ) = 0 

a - x' 


<‘a  - V ■ I5o  C1  - + \ + **2>  - * 10-4 


Combined  effect  on  market  share  z of  both  approximations 

The  correction  terms  for  the  combined  effect  on  z of  both  model- 
ing approximations  are: 


A z 

■ (z  - z ')  + 
a x ' 

(z  - z ) 
x m 

V 

V V 

V V 

A z 

“ (z  - z ) + 

a x 

(z  - z ) 
x m' 

160 

, these  correction  terms 

are 

A z 

- -.0011 

+ 0 

V 

-4 

A z 

- 3.30  x 10 

+ 1.86  x 

10" 

-.0011 


= 5.16  x 10 
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Then,  the  estimated  authentic  mean  and  variance  of  z for  p ■ 160 
are: 

z = z + az  = .846 

am 

V V V _A 

z = z + Az  = 43.61  x 10 

a m 

We  see  that,  for  p = 160  , the  analysis  overstates  the  mean  of  z 

by  about  0.1%  and  understates  the  variance  of  z by  about  12%  and 
the  sta.  .rd  deviation  by  about  6%. 


Effect  on  market  size  n of  the  Distribution  Approximation  on  the 
economic  indicator  f : 

From  the  Demand  Model,  we  have: 


, f ( 225 
n = Jjf  e 

Taking  derivatives  with  respect  to  f : 


n'(f) 


he 


-Ilfs 


n"(f>  = 0 

For  p » 160  , x = 160  so: 

n' (f)  = .246  for  all  f 

Noting  that  n"(f)  =0  , we  have  the  correction  terms  for  the  effect 

on  n of  the  Distribution  Approximation  from  (2.11)  and  (2.12)  : 
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(a  “ n ) « n'(£  ) Af 

x in  m 


V V 9 v 

(nx  • V ■ n VJ  ' Af 


Substituting  from  Table  3.3  the  values  of  Af  and  Af  : 


■ (-0.014)  n’(fm) 


(n  - n ) - (0.005)  n'(f)' 

X IQ  tu 


. - 2 


For  p ■ 160  , t'neje  correction  terms  are: 


(n  - a ) - -.0035 

X IQ 


(n  - n ) **  3.24  x 10^ 

X IQ 


Effect  on  market  size  n of  the  Dependence  Approximation  in  the 
Demand  Model: 

The  Demand  Model  claims  that  only  the  lower  celling  price  and  the 

general  economic  conditions  affect  the  size  of  the  total  market.  As 

with  the  Market  Split  Model,  the  client  believes  that  other  factors  may 

also  influence  the  total  market  size  and  that,  therefore,  the  model 

understates  his  uncertainty  about  the  market. 

We  characterize  this  Dependence  Approximation  by  the  random 

variable  e : 
n 
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The  client  assesses  his  residual  uncertainty  as: 


<ep|x,f,4>  = 0 

v<ejx,f,4>  = (1.5  x 10_5)(x  - 120)  f2 


Again,  since  the  conditional  mean  is  zero,  e * 0 and 

n 

v 

cov(nx,en)  * 0.  Then,  calculating  e^  : 


v<en|x,M){f  I*} 


(1.5  x 10"5)  (x  - 120)  | f2{f|<0 


(1.5  x 10"5)  (x  - 120)  f 2 


» (1.5  x 10"5)  (x  - 120)  (f  + f 2) 

3 3 


Substituting  from  Table  3.3  the  values  of  f and  f : 

3 3 

e = (1.55  x 10"5)  (x  - 120) 


The  correction  terms  for  the  effect  on  n of  the  Dependence  Approxi- 
mation from  (2.14)  and  (2.15)  are: 


(n  - n ) = e 
ax  n 


v v v v 

(n  - n ) = e + 2cov(n  ,e  ) = e 
a x n x’  n n 


For  p ■ 160  , these  correction  terms  are: 


(n  - n ) = 0 
a x 

v V — A 

(n.  - n ) = 6.19  x 10 

3 X 


> - VT 1 1 1'Js1  ¥ W 'M j.'.» . JU.'.IU ^ l1* W 1 I'  mi hiwuiiii (II JU1^1].|1, j^vi^uijpf  jip  %-u ! ■ 


Combined  effect  on  market  size  n of  both  approximations: 

The  correction  tenu  for  the  combined  effect  on  n of  both  model- 
ing approximations  are: 


An 

v 

An 

For  p * 160 


(n  - n)  + (n  - n) 
a x7  x m 


v v v v 

(n  - n ) + (n  - n) 
a x x m7 


these  correction  terms  are: 


An  » -.0035  + 0 = -.0035 

An  = 3.24  x 10~4  + 6.19  x 10-4  - 9.43  x 10 


The  it  deled  mean  and  variance  of  n for  p = 160  are: 

n - .250 
in 

n - 11.95  x lO-4 
m 

So,  the  estimated  authentic  mean  and  variance  for  p = 160  are: 

n - n + An  ■ .246 

a m 

V V V _A 

n * n + An  * 21.38  x 10 

a m 

We  see  that,  for  p ■ 160,  the  analysis  overstates  the  mean  of  n by 
about  1.42  and  understates  the  variance  of  n by  about  44%  and  the 
standard  deviation  by  about  25%. 
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Cumulative  effect  of  the  modeling  approximations  on  total  sales  q : 
Now,  considering  the  output  variables  of  the  Market  Split  and 
Demand  Models,  market  share  z and  total  market  size  n , as  input 
variables  to  the  next  submodel,  we  can  determine  the  cumulative  effect 
on  total  sales  q of  the  various  modeling  approximations.  These 

v 

approximations  are  characterized  by  the  correction  terms  A z , A z , 
- v 

An  and  An  that  we  have  just  computed. 

From  the  model,  we  have: 

q = z n 

This  is  an  exact  relationship,  so  there  is  no  Dependence  Approximation 
for  this  submodel;  we  need  worry  about  only  the  Distribution  Approxi- 
mation on  input  variables  z and  n.  Taking  partial  derivatives: 

q (z,n)  - n q (z,n)  - 0 

z zz 


qn(z»n)  - 2 qnn(z,n)  “ 0 

qzn(z-”>  ‘ 1 

For  p ■ 160  : 

’A.V  ■ -250 

VVV  = -847 

Noting  that  qzz  and  qnn  are  zero,  we  have  the  correction  terms  for  the 
effect  on  q of  the  Distribution  Approximation  from  (2.29)  and  (2.30) 
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(q  " q_)  “ q (z  ,n  )Az  + q (z  ,n  )An 
Mx  m z m m n m m 


+ ^„„(z  »n„ )AzAn 

zn  m m 


v v 2 ' _ _ 2 v 

(q  “ q_)  * q (z  ,n  ) Az  + q (z  ,n  ) An 
Tn  z m m n in  m 


_ _ 2VV  v v v _2  v _2 
+ q (z  ,n  ) (z  n -z  n +z  An  +n  Az  ) 
zn  m m aamma  a 


+ 2q  (z  ,n  )q  (z  ,n  ) z An 
Mz  m m Mzn  m m a 


+ 2q  (z  ,n  )q  (z  ,n  ) n Az 
n m m zn  mm  a 


Since  there  is  no  Dependence  Approximation: 


a5  - (5a  - 5„)  + <ix  - qa)  - Cqx  - V 


V V 


v v 


v v 


Aq  - (qa  - qx)  + (qx  - q„>  - (q„  - ij 


For  p ■ 160  , these  correction  terras  are : 


Aq  - -.0032 


Aq  - 7.02  x 10~4 
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The  modeled  mean  ana  variance  of  q for 


P ■ 160  are: 


Sn  " *212 


qm  - 11.01  x 10' 


So,  the  estimated  authentic  mean  and  variance  o2  total  sales  q 


for  p «=  160  are: 


*a  ” \ + ■ -209 


\ + Aq  - 18.03  x 10" 


We  see  that,  for  P * 160  , the  analysis  overstates  the  mean  of 

q by  about  1.5%  and  understates  the  variance  of  q by  about 
39%  and  the  standard  deviation  by  about  22%. 


Cumulative  effect  of  modeling  approximations  on  profit  s : 

In  the  model,  profit  s is  given  as  the  following  function  of 
selling  price  p , total  sales  q , fixed  costs  t , other  variable 
costs  b , raw  material  cost  h and  material  usage  factor  w : 

s=pq-t-q(b+  . 85hw  + .15hw  e"qW) 

= q (P  - b - . 85hw  - . 1 5hw  eKqW)  - . 


where  K = -ln(81) 

This  is  an  exact  relationship,  so  there  is  no  Dependence  Approximate 
Taking  partial  derivatives: 


-wq ( . 85  + .15  e qW) 


sw  * -hw(.85  + .15(1  + Kq)  eqw) 

3q  = P - b - hw(,85  + .15(1  + Kqw)  eKqW) 


at  = _1 


8bb  " 0 


Shh  * 0 


-.15  Khq  (1  + Kq)  e 


sqq  = -.15  Khw2  (2  + Kqw)  eKpW 


Stt  = 0 


Sbh  = 0 


s.  ■ 0 

bw 


“bq  - -1 


s(.)  - 

Then,  for  p - 160  , the  partial  derivatives  evaluated  at  the 

modeled  mean  are: 


8b(.) 

* 

-.212 

8h(.) 

00 

co 

CM 

• 

i 

8w(’) 

-7.072 

8q(* ) 

- 

78.353 

8ww^’  ^ 

*3 

.025 

03 

• 

v-/ 

- 

10.500 

/-% 

s-/ 

i= 

CO 

- 

i 

• 

00 

8h,(>) 

“ 

-1.053 

8wq'-> 

S3 

-31.509 
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Noting  that  Ay  , Aw  and  At  are  all  zero,  we  have  the  correction 
terms  for  the  effect  on  s of  the  Distribution  Approximation  from 
(2.27)  and  (2.28)  : 


s (.)Ab  + s (.)Ah  + s (.)Aq 
D n q 

+ h s (.)Aw  + Is  s ( . ) (Aq  +Aq  2) 
ww  qq 

+ 8^  (.)AbAq  + s^(.)AhAq 

« v o * 2 v 

8.  (.)  Ab  + s.  (,)  Ah  + s (.)  Aw 
d n w 

2 v 2 v , 2 v2 

+ 8 (.)  Aq  + s (.)  At  + **8  (•)  Aw 

q t ww 

. 2 w 2 v 9 v -2 

+ h 8 (•)  (qfl  -q  ^ 2q  Aqz) 

qq  & in  & 

oVV  vv  v ov_o 

+ 8bq<0  (baqa  ' Vm  + Vq  * qaAb  ) 

2 v v vv  v _2  V _ 2 
+ s.  (.)  (h  q -hq  + hAq  + q Ah  ) 
hq  a a mm  a a 

2 v _o  v 2 v _2  v 

+ shw(.)  (hfl  + Ah  )Aw  + swq(.)  (qfl  + Aq  )Aw 

+ 28bq(-)8hq(-)qaAbAb  + 8qq(* )8bq(' )qaAqAb 

+ 8qq(-)8hq(')qaAqAb  + 28b(*)8bq(* )baAq 

4 28h(,)8hq(,)haAq  + 28w(*)8hw('>^Ab 

+ 28w(.)s£iw(.)AwAq  + 2sq(.)sq(?(.)qaAq 

+ 2s  ( • ) 8,  ( . )q  Ab  + 2s  (.)s  (.)q  Ah 

q bq  a q qh  a 
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because  there  is  no  Dependence  Approximation: 


As  - (s  - s ) + (s  - s ) - (s  - s ) 

ax  xm  xm 


As  » (s  - s ) + (s  - s ) - (s  - s ) 

“ ' a v'  V flj'  v m 


m 


For  p ° 160  , these  correction  terms  are: 


As  * -.472 


As  - 5.797 


The  modeled  mean  and  variance  of  s for  p ■ 160  are: 


s = 7.014 

m 

s = 11.317 

m 

So,  the  estimated  authentic  mean  and  variance  of  profit  s for 
p = 160  are: 


s 

- 8 

+ As  “ 

6.542 

a 

m 

V 

V 

V 

s 

- s 

■1-  As  -- 

17.114 

a 

m 

We  see  that, 

for  p 

= 160  , 

the  analysis  overstates  the  expected 

profit  by  about  7.2%  and  understates  the  variance  by  about  34%  and 
the  standard  deviation  by  about  19%. 


Cumulative  effect  of  the  modeling  approximations  on  the  certain  equ‘ .va- 
lent : 

By  treating  the  utility  function  as  a submodel,  we  can  determine 
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the  emulative  effect  of  all  the  modeling  approximations  on  the  certain 
equivalent. 

The  utility  function  is ; 
u(8)  - 1 - e ^ 8 

Taking  derivatives: 

u’(s)  - .1  e~‘l  * u"(s)  - -.01  e"'1  8 

For  p ■ 160  , these  derivatives  evaluated  at  s are: 

in 

u*(s  ) - .0496 
m 

u"(s  ) - -4.96  x 10~3 
m 

The  correction  term  for  the  effect  on  expected  utility  of  the  Distri- 
bution Approximation  from  (2.11)  is: 

Au  ■ u' (s  )Aa  + u’’(a  ) (As  +As^) 

m m 

For  p ■ 160  , this  correction  term  is : 

Au  - -.0383 

The  modeled  expected  utility  for  p ■ 160  is: 

u - .4764 

m 

So  the  estimated  authentic  expected  utility  for  p - 160  is: 

u - u + Au  - .4381 

a m 


90 


This  corresponds  to  a certain  equivalent  of : 

CE  « u"1  (uj  - 5.76 

The  modeled  certain  equivalent  for  p ■ 160  is: 

CE  - 6.47 

m 

We  see  that,  for  p ■ 160  , the  modeling  approximations  have  caused 

the  analysis  to  overstate  the  certain  equivalent  by  about  $710,000 
or  about  12%. 

By  repeating  this  procedure  for  different  values  of  p , we 
find  the  following: 

Modeled  Est.  authentic  Overstatement 


p 

CE 

m 

CE 

a 

by  analysis 

150 

5.87 

5.20 

13% 

160 

6.47 

5.76 

12% 

170 

6.40 

5.67 

13% 

After  taking  into  account  the  effect  of  the  modeling  approximations,  we 
see  that  the  optimal  decision  is  still  p ■ 160. 

Effect  of  the  modeling  approximations  on  the  value  of  clairvoyance: 

If  we  apply  the  methodology  to  the  probabilistic  sensitivity 
analysis,  we  can  estimate  the  effect  of  the  modeling  approximations  on 
the  value  of  clairvoyance.  We  find  as  a result  that  the  analysis 
slfghtly  overstates  the  value  of  clairvoyance: 
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Value  of  Clairvoyance 


Variable 

Modeled  Est.  authentic 

Competitor's  price: 

y 

$378,000 

$327,000 

Economic  indicator: 

f 

0 

0 

Raw  material  cost: 

h 

$ 46,000 

$ 43,000 

Other  var.  costs: 

b 

$ 31,000 

$ 25,000 

Summary 

We  see  that  the  methodology  presented  in  Chapter  2 can  be  used  to 
determine  the  effect  on  the  profit  lottery  and  on  the  certain  equivalent 
of  the  modeling  approximations  made  in  the  analysis.  The  methodology 
thus  provides  us  with  an  indication  of  the  adeqjacy  of  the  modeling  in  a 
decision  analysis. 

In  this  example,  the  modeling  appears  to  be  quite  adequate,  for  the 
methodology  indicates  that  the  modeling  approximations  do  not  affect  the 
optimal  decision. 


CHAPTER  4 


STOCHASTIC  MODELS 


t 


4.1  Introduction 

In  the  preceding  chapters,  we  have  dealt  solely  with  deter- 
ministic models  - those  models  that  yield  a unique  value  of 
the  outcome  variable  for  fixed  values  of  the  state  and  decision 
variables.  We  can  now  look  at  stcchastic  models  - those  models 
that  yield  a probability  dts;,,hution  on  the  outcome  variable  for 
fixed  values  of  ti  e state  ai.u  dacision  variables. 

We  use  stochastic  models  to  represent  situations  in  which 
we  are  uncertain  about  the  dependence  of  the  outcome  variable  on 
the  state  variables  but  cannot  conveniently  attribute  this  uncer- 
tainty to  additional  state  variables.  Indeed,  in  these  situations, 
we  often  view  the  dependence  relationship  as  being  inherently 
uncertain.  Generally,  we  possess  some  data  (e.g.,  long-run 

frequencies)  about  the  uncertain  relationship  to  help  us  define  the 
stochastic  model. 

For  example,  we  use  the  Bernoulli  stochastic  model  to  express 
our  uncertainty  about  the  number  of  heads  obtained  in  n tosses  of 
a coin  because  we  find  it  easier  to  think  of  the  outcome  of  a coin 
toss  as  being  inherently  uncertain  than  to  think  of  it  as  being 
dependent  on  such  extraneous  factors  as  rotational  forces  and 
surface  characteristics. 
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In  this  chapter,  we  shall  extend  the  methodology  developed 
in  Chaptei  2 for  deterministic  models  to  stochastic  models.  The 
methodology  is  useful  only  for  relatively  simple  stochastic 
models,  such  as  Bernoulli  or  small  Markov  models. 

4.2  Extension  of  the  Methodology 

In  a decision  analysis,  we  ideally  want  the  authentic  profit 
lottery: 

{v|d,4}  - | {v|d,s,4}{sj4} 

5. 

Suppose,  however,  that  we  use  a stochastic  model  to  compute  the 
profit  xottery,  as  follows: 

(v)tl,4}m  - l {v|d,s.,4}  p(s) 

£ 

where  the  stochastic  model  is  specified  by  the  conditional  proba- 
bility distribution  {v|d,8,4)m.  We  see  that  the  discr»pat.cy 
between  the  modeled  and  authentic  profit  lotteries  is  due  to  the 
two  modeling  approximations: 

1.  The  approximation  of  the  probability  distribution  on  the 
state  variable  vector  s by  the  discrete  probability  function  p 
(Distribution  Approximation) : 

(4.1)  p(s)  = {s>} 


94 


2.  The  approximation  of  the  conditional  probabllit> 
distribution  on  v given  d and  s^  by  the  stochastic  model 
(Dependence  Approximation) : 


(4.2) 


{vld,s,i} 


{v|d,s,i} 


We  want  to  quantitatively  characterize  these  modeling  approxi- 
mations and  determine  their  effect  on  the  profit  lottery. 

More  specifically,  letting  vq  denote  the  authentic  out- 
come variable  and  vy  denote  the  outcome  variable  of  the 
stochastic  model,  we  want  to  determine  the  correction  terms: 


(va  - vy) 


V V V 

Av  „ (va  - vy) 


Let  us  first  define: 


g (d , s)  a <vy|d,s,4> 


f(d,s)  = V<vy!d,s,4> 


the  conditional  mean  and 
variance  of  the 
stochastic  model 


Then,  we  can  write: 


(4.3) 


g(d,s)  + £ 


' 


where  e is  a random  variable  whose  conditional  mean  and 


variance  are: 


<e|d,s,4> 


v<e|d,s,  i) 


(vyld»£»4)  “ g(d»s)  - 0 


V<vyld»s,4> 


f(d,s) 


Then,  using  g(d„s)  as  a reference  deterministic  model, 
we  can  write  the  Dependence  Approximation  in  (4.2)  a a two 
separate  approximations : 


Kjd.s,*}  2 5(v  - g(d,s))  = (v  |d,s,A} 
y m a — 

We  characterize  the  first  of  these  approximations  by  the  random 
variable  e as  in  (4.3)  above  .nd  the  second,  as  in  Chapter 
2,  by  the  random  variable  e : 


(4.4)  v = g(d,s)  + e 


where  the  conditional  mean  <c|j,£,4>  and  variance  v<e|d,£,A> 
are  directly  assessed. 

Next,  employing  the  following  notation  from  Chapter  2: 


n - the  state  variable  vector  whose  distribution  is  p(s) 
x ■ the  state  variable  vector  whose  distribution  is  {sli} 
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we  characterize  the  Distribution  Approximation  in  (4.1)  by 
the  correction  terms: 


x - n 


v V V 

As  * x - n 


We  use  four  different  subscripts  for  the  outcome  variable 
(see  Figure  4.1): 

v ■ the  outcome  variable  of  the  stochastic  model 

y 


v * the  outcome  variable  of  the  reference 
m 

deterministic  model 

v ■ the  outcome  variable  corrected  for  the 
x 

Distribution  Approximation 


v = the  authentic  outcome  variable 
a 


Then,  we  can  write: 

v - g(d,n) 
m 


V * g(d,x) 
x 


(4.5)  vy  - g(4;n)  + e = vm  + F- 


(4.6) 


v - g(d,x)  + e - -r  + e 


Taking  the  mean  and  variance  of  (4.5),  we  have: 
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Figure  4.1. — Definition  of  Subscripted  Omrrm»> 


Vy  " Outcome  variable  of  stochastic  model 


{vy'd,4}  - l {v|d,s,4>  p(s) 


Correction  for: 

{v|d,s,4}ffl  = <s(v  _ g(d>s)) 


Vm  “ variable  of  reference 

I deterministic  model 


{vuJd’4}  - l 6(v  - g(d,s))  p(s) 


Correction  for: 


vx  - Outcome  variable  corrected  for 
Distribution  Approximation 

{vjd.4}  - f 6(v  - g(d,s)){s|4} 

I • a 


Correction  for: 

6Cv-g(d,s))  * {v|d,s,4} 


va  * Authentic  outcome  variable 

{vjd,4}  - [ {v|d,s,4}{s|4} 


«T», 


V V V 

v ■ v + e + 2cov(v  ,e) 

y m m 

But,  because  <e|d,s,4>  - 0 , e - 0 and  cov(v  ,e)  - 0. 

m 

So, 

(4.7)  v - v 

y m 

V V V 

(4.8)  v - v + e 

y m 

where: 

£ ■ l V<e|d,s,4>  p(s) 


(4.9) 


- I f(d,s)  p(s) 
8 


Next,  taking  the  mean  and  variance  of  (4.6),  we  have: 

(4.10)  v 


(4.11) 


v + e 
x 


V V 

v + e + 2cov(v  ,e) 

^ X 
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where  (see  Appendix  A. 4)  : 


(4.12) 

e “ 

• 

<e|d,s,4){sj4} 

(4.13) 

V 

e ■ 

J 

(v<e|d,£,4)  + <e|d,£,4)2){sj4}  - e2 

s_ 

(4.14) 

cov(v^,e] 

1 - J g(d,8)<e|d,8_,4){sj.6}  - v^e 

5.  X 

Then,  subtracting  (4.7)  from  (4.9)  and  (4.8)  from 
(4.10),  we  have: 

(4.15)  Av  ■ (v  - v ) * (v  - v ) + e 

a y x m 

v v v v v v v 

(4.16)  Av  - (vfl  - Vy)  ■ (v^  - v^)  + (e  - e)  + 2cov(vx,e) 

v _ v 

where  e , e , e and  cov(vx,e)  are  calculated  as  in  (4.9) 

and  (4.12)  through  (4.14)  above  and  where  (v  - v ) and 

x m 

v v 

(v  - v ) are  the  correction  terms  for  the  effect  of  the  Distri- 
x m 

bution  Approximation,  given  by  (2.11)  and  (2.12)  for  the  single 
state  variable  case  and  by  (2.27)  and  (2.28)  for  the  multi- 
variate case  The  '■-.quations  (4.15)  and  (4.16)  above  yield 
the  desired  correction  terms  for  the  outcome  variable. 

As  a special  case,  suppose  that  the  stochastic  model  fully 
expresses  our  uncertainty  about  the  dependence  relationship: 


^Vald»£»-6^ 
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r 


Then, 


<v  |d,s,4> 

a 


<vy|d,s,4> 


so, 


[ 


(e | d,£,4) 

and,  therefore. 


<v  Jd,s,4>  - g(d  ,s) 


0 and  covCv^je) 


Also, 


V(e|d,£,4)  - ^Vyld.^,4)  - f(d,s) 


and,  therefore,  from  (A. 13), 

(A. 17)  e - [ f(d,s){s|A} 

* R 


Then,  from  (A. 9)  and  (A. 17), 


v v 
( e - e)  - 


J 


f(d,s){8|4}  - l f (d,s)  p(jO 


— ®. 

From  (2.33),  we  see  that  this  can  be  estimated  as: 


(e  - e)  = l f (d,n)Ai  +hl  f..(d,n)As, 
1 1 11  1 


(4.18) 


*5  I I f,  .(d,n)As  As 
1 i 1J  1 J 
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So,  for  this  special  case,  the  correction  terns  in  (4.15) 
and  (4.16)  become: 

(4.19)  Av  ■ (v  - v ) 

x m 

v v v v v 

(4.20)  Av  - (v  - v)  + (e-e) 

x m 

4.3  Example:  A Game 

Consider  the  following  game:  We  draw  one  ball  from  an  urn 

containing  25  balls  that  are  numbered  consecutively  from  1 to 
25.  Letting  n be  the  number  on  the  drawn  ball,  we  then  flip 
a thumbtack  n times.  Finally,  letting  r be  the  number  of  "heads" 
(point  up)  obtained,  we  receive  a prize  of  v dollars,  where  v - r2. 

Should  we  pay  $25  to  play  this  game?  Our  utility  function  on 
dollars  is  exponential  with  risk  aversion  coefficient  y » .04  : 

«(v)  - 1 - .'•04v 

We  use  the  Bernoulli  stochastic  model  to  analyze  this  game. 

(See  Figure  4.2)  The  state  variables  are  n , the  number  of  flips, 
and  p , the  long-run  frequency  of  "heads"  in  many  flips  of  the 
thumbtack.  Given  fixed  values  of  these  state  variables,  the 
Bernoulli  model  yields  the  following  probability  distribution  on 
the  number  of  "heads"  r : 

{r|n,p,4}m  - j"](p)r(l  - p)n"r 
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..  .-a.  -t  -i  - -1-  - — •- 


n - number  on  drawn  ball 
P - long-run  frequency  of  heads 
r ■ number  of  heads  in  n tosses 
v ■ monetary  payoff 
u « utility 
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For  the  purposes  of  analysis,  we  use  the  following  distributions 
on  the  state  variables: 


{n|4} 


m 


2^  n - 1,2, ...25 
0 otherwise 


{p|4> 


m 


<5<P  “ P0)  where  pQ  - .5 


Using  a computer,  we  calculate  the  resulting  probability 
distribution  on  r : 


{r|4} 


m 


l £ {r|n,p,4}m{n|4}m  {p|4)m 
n p 


Next,  we  calculate  the  corresponding  distribution  on  v : 


H4}m  " l {v|r,4}{r|4}m 
r 

- I 6(v  - r2)  {r |-6 > 


Finally,  we  compute  the  expected  utility  and  the  certain  equivalent 
of  the  game: 


u ■ 

m 


j>  (v 1 4 } u(v) 


m 


CE  - u_1(u  ) 

m m 
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We  find  that: 


v 


r 

m 

- 6.50 

r ■ 

m 

16.25 

V 

n 

- 58.50 

V 

v ■ 

m 

3586.05 

u 

m 

- 0.6469 

CE 

- 26.03 

m 


The  analysis  indicates  that  we  should  be  willing  to  pay  up  to 
$26.03  to  play  this  game. 

We  can  now  ask:  What  modeling  approximations  have  we  made 

in  this  analysis  and  what  effect  do  they  h.:ve  on  its  results? 

In  this  case,  there  is  just  one  approximation  — the  use  of  a 
fixed  value  p^  • .5  for  the  long-run  fraction  of  "heads'* 
instead  of  a probability  distribution  on  p.  Suppose  that  we 
assess  the  following  mean  and  variance  for  p : 

P * P0  - -5 

v 

p = .01 

Then,  we  characterize  the  approximation  by  the  correction  terms 
Ap  * 0 

Ap  - .01 
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model  the  conditional 


We  define  as  the  reference  deterministic 
mean  of  the  Bernoulli  .1,  ie-.t‘ng: 

g(n.p)  - (r |n,p,a)  - n p 

f(n,p)  ■ (r|n,p,4)  - n p (1  - p) 

Taking  partial  derivatives  of  g(n,p)  : 


Then,  noting  that  An 
(2.27)  and  (2.28)  : 


nn 


PP 


np 


v 

An  and 


■ 0 

- 0 
- 1 


Ap 


are  all  zero,  we  have  from 


0 

8 (n,p)2Ap  + g np(n,P)2(p  n - ^ n) 
r ur  a m 

-2  V 7 V 
n Ap  + nAp 


-2  v v 
(n  + n)Ap 


Because  the  Bernoulli  model  fully  expresses  our  uncertainty 
•bout  the  number  of  ■■head."  given  n and  p , v,  can  use  the 
results  of  Section  4.2  for  the  special  case.  Taking  partial 
derivatives  of  f(n,p)  with  respect  to  p : 


fp  “ n - 2np 


f __  “ -2n 

PP 


Then,  from  (4.18),  we  have: 


v v 
(e  - c) 


^ppCr.,p)Ap 


— V 

-n/ip 


So,  from  (4.19)  and  (4.20),  we  have: 


a; 

■ (r  - r ) 
x m 

- 

0 

V 

Ar 

V V v 

■ <r*  * r„>  + <e 

V 

“ e)  ■ 

(n‘2  + 

Substituting,  we  find  that: 
Ar  * 0 

*a  " 

6.50 

V 

Ar 

- 2.08 

V 

r ■ 

a 

18.33 

- v 
n)Ap 


So,  the  fixation  of  p has  no  effect  on  the  expected  value  of  r 
but  causes  the  analysis  to  understate  its  variance  by  about  11Z. 

We  can  now  determine  the  effect  of  the  approximation  on  our 
winnings  v.  Taking  derivatives  of  v ■ r^  : 


2r 

2 
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Then,  noting  that  Ar  - 0 , we  have  from  (2.11)  and  (2.12): 


I 


Av  ■ *jv"(r  )Ar 

tn 

v 

- Ar 

= - 2 v _ •>  v v v 

Av  - v’  (r  ) Ar  + %v"(r  ) Ar(r  f r ) 

n mam 

_ 9 v v v v 

- Ar  Ar  + 2Ar(2r  +Ar) 

m m 

Substituting,  we  find  that: 


Av  - 2.08 

V 

a 

- 60.58 

V 

V 

Av  - 495.37 

V 

a 

- 4081.42 

We  see  that  the  fixation  of  p causes  the  analysis  to  understate 
both  the  expected  value  and  the  variance  of  v by  about  3. / and 
12Z  , respectively. 

Finally,  we  can  determine  the  effect  of  the  approximation  on 
the  expected  utility  and  the  curtain  equivalent.  Taking  derivatives 
of  the  utility  function: 
u(v)  - 1 - 

u'  - .04e“-04v 

u"  - -.0016  c ■•04v 
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Then,  from  (2.26),  we  have: 


Au  - u' (v  )Av  + *su"(v  ) (Av  +Av  ) 
in  n 

■ .04  e *^Vm  (Av  - .02(Av  +Av^)) 

Substituting,  we  find  that: 

Au  - -0.0305  ua  - 0.6164 

a 

And: 

CE  - -25  ln(l  - u ) - $23.96 

a a 

So,  we  see  that  the  fixation  of  p causes  the  analysis  to  over- 
state the  certain  equivalent  by  about  8.6%  and  that  we  should 
not  be  willing  to  pay  $25  to  play  the  gamt. 

It  is  interesting  to  see  how  the  misstatement  of  the  certain 
equivalent  caused  by  the  fixation  of  p varies  with  the  degree  of 
our  risk  aversion.  Figure  4.3  showp  the  percentage  misstatement 
of  the  certain  equivalent  as  a function  of  our  risk  aversion 
coefficient  Y for  the  game  played  with  5 , 15  and  25  balls 
in  the  urn.  Note  that  the  misstatement  of  the  certain  equivalent 
may  be  either  positive  or  negative,  depending  on  our  risk  preference, 
and  that  the  degree  of  misstatement  does  not  necessarily  vary  mono- 
tonically  with  the  risk  aversion  coefficient. 
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CONCLUSION 


As  we  have  seen,  the  role  of  modeling  in  decision  analysis 
is  to  help  us  obtain  the  authentic  profit  lottery  by  simplifying 
our  assessment  task.  We  think  of  the  models  in  decision  analysis, 
then,  as  subjective  expressions  of  our  uncertain  understanding  of 
the  world  rather  than  as  objective  descriptions  of  reality. 

We  can  new  answer  the  questions  about  modeling  that  we 
raised  in  the  Introduction. 

First,  what  is  "goodness"  in  a model  and  can  we  quantify  it? 
Becausa  we  tnink  of  a model  as  an  approximate  expression  of  our 
uncertainty  and  not  as  a description  of  the  real-world,  an  index 
of  realism  would  not  provide  a meaningful  measure  of  goodness. 
Rather,  the  measure  of  goodness  that  we  want  is  how  well  the  model 
represents  our  uncertainty  about  the  dependence  relationship  between 
outcome  and  state  variables.  We  can  quantify  this  measure  of  good- 
ness with  the  random  variable  e : 

v - g(d,ai)  + e 

Then,  we  can  use  the  methodology  of  Chapter  2 to  show  how  the  good- 
ness of  the  model  affects  the  results  of  the  decision  analysis  by 

determining  its  contribution  to  the  correction  terms  on  the  profit 

v 

lottery  £v  and  Av. 


Ill 


^ - 


Next,  how  should  we  choose  among  alternative  models?  Using 
the  measure  of  goodness  defined  above,  we  should  choose  among 
competing  models  according  to  how  ’'close"  they  get  us  to  the 
authentic  profit  lottery.  When  our  state  of  information  about 
the  decision  environment  includes  a large  amount  of  data,  we  can 
use  a probability-space  of  models  to  help  us  choose  among  several 
alternative  model  forms,  being  careful  that  we  meaningfully 
interpret  the  probability  assigned  to  each  model  in  the  space. 

(See  Section  1.6.) 

Finally,  how  should  we  decide  when  to  do  more  modeling?  The 
reason  that  we  use  modeling  in  decision  analysis  at  all  is  that  we 
judge  the  directly  assessed  profit  lottery  to  be  so  "far"  from  the 
authentic  profit  lottery  that  it  is  worthwhile  spending  the  time 
and  effort  to  reduce  our  secondary  uncertainty.  Likewise,  our 
decision  to  do  more  modeling  should  depend  on  whether  or  not  it 
is  worth  the  additional  effort  to  £~t  even  "closer"  to  the  authentic 
profit  lottery. 

We  use  the  concept  of  the  value  of  perfect  information  to  help 
us  decide  when  to  collect  more  data.  The  analogous  concept  that  we 
need  to  help  us  decide  when  to  do  more  modeling  is  not  the  value  of 
perfect  modeling  but  rather  the  value  of  obtaining  the  authentic 

lottery.  This  value  is  the  most  we  should  pay  to  completely 
resolve  our  secondary  uncertainty  through  further  modeling.  Unfor- 
tunately, it  is  exceedingly  difficult  to  exercise  this  concept,  since 
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to  do  so,  we  must  assess  the  joint  distribution  on  the  authentic 
certain  equivalents  for  all  decision  alternatives. 

We  can,  however,  get  an  indication  of  the  value  of  additional 
modeling  by  using  the  methodology  of  Chapter  2 to  estimate  the 
difference  between  the  modeled  and  authentic  profit  lotteries.  The 
greater  the  difference,  the  greater  the  presumed  value  of  additional 
modeling.  Furthermore,  we  can  use  the  methodology  to  indicate  where 
more  modeling  would  be  most  effective  by  showing  the  contribution  of 
each  part  of  the  model  to  the  discrepancy  in  the  profit  lotteries. 


appendices 


Appendix  A.l:  Effect  of  the  Distribution  Approximation 

This  appendix  supplements  Section  2.2,  pp.  39-40. 

For  notational  ease,  let: 

g(.)  ■ g(d,n) 

3 - 3 

8 “ (s  - s)  , third  central  moment 

8 (s  - 8)  , fourth  central  moment 

We  characterize  the  Distribution  Approximation  by: 

As  = x - n 

v v v 

As  = x - n 

We  expand  vm  - g(d,n)  about  n : 

(a.l)  vm  * g(.)  + g'(.)  (n  - n)  + ?ig"(.)  (n  - n)2 

Taking  the  mean: 

(a. 2)  vffl  - g(.)  + isg"(.)  n 

Squaring  (a.l)  and  taking  the  mean: 

vm2  = 8(.)2  + g'(.)2  (n  - n)2  + *sg"(.)2  (n  - n) 

+ 2g(.)g'(.)  (n  - n)  + g(.)g"(.)  (n  - n)2 
+ g'(.)g"(.)  (n  - n)3 
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(a.3)  vm2  = g(.)2  + g'(,)2n  +JSg"(02  n 


+ g(.)g’'(.)n  + g'(.)g"(.)  n 


Squaring  (a. 2)  : 

(a.4)  vffi2  = g(.)2  +Hg"(.)2  n2  + g(.)g"(.)  n 
Subtracting  (a.4)  from  (a.3)  : 

(a,5)  vm  = 8’(*)2  n + Jsg"(.)2  (n  - n2)  + g'(.)g"(.)  n 

To  simplify  this  expression,  we  assume  the  following  (see  Appendix 
A.  2)  : 

3 

n ■ 0 


n " 3 n 


Then  (a. 5)  becomes: 


(a. 6)  vm  a g*(.)2  I + >sg"(.)2  n2 


Similarly,  we  expand  - g(d,x)  about  ii  : 

(a*7)  vx  " 8(.)  + g’(.)  (x  - n)  + WC.)  (x  - n)2 


Taking  the  mean: 


(a. 8)  vx  « g(.)  + g'(.)  (x  - n)  + Jsg"(.)  (x  - n)2 


' 


Squaring  (a, 7)  and  taking  the  mean: 

v 2 = g(02  + g'(-)2  (x  - n)2  + Hg"(*)2  (x  - n) 


(a. 9) 


-v  2 


+ 2g(.)g'(.)  (x  - n)  -*•  g(.)g"(.)  (x  - n) 
+ g' (.)g"(.)  (x  - n)3 


-v  2 


-.4 


v 2 = g(.)  + g'(.)  (x  - n)  + W(.)  (x  - n) 


-v  2 


+ 2g(.)g'(0  (x  - n)  + g(.)g"(.)  (x  - n) 


+ g'(.)g"(.)  (x  - n)3 


Squaring  (a. 8)  : 


(a. 10) 


v 2 = g(.)2  + g'  (.)2  (x  - n)2  + ^<g"(.)Z  (x  - n) 


2 


-.2 


+ 2g(.)g'(.)  (x  - n)  + g(.)g"(*)  (x  - n) 


-.2 


+ g'  (Og"(.)  (x  - n)  (x  - n) 


Subtracting  (a. 10)  from  (a.9> 


(a. 11) 


g-  (.)*!  (X  - n)2  - (x  - n)2 


!] 


+ g'( 


L 

.)g"(*)  (x  - 


] 


+ *sg"(.)2  |jCx  - n)^  - (x  - n)2 

n)3  - (x  - n)  (x  - n)2^ 
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Note  that; 


(x  ~ n) 


As 


(x  - n)2 


(x  - il)3 


(x  - n)^ 


( (x  - x)  4*  (x  - n))  ■ x + As  ^ 


((x  - x)  + (x  - n)) 


rss3 


3 V _ 
x + 3xAs+As J 


((x  - x)  + (x  - n))  « x + 4xAs+6xAs2+  As 


v _2 


Substituting  into  (a. 8)  and  (a. 11)  : 

(a. 12)  v;£  « 8(0  + g'(.)As  + >jg"(.)  (x+As 2 ) 


'a. 13) 


„ 2 v 


= g'  (•)"  x + %g"(.)*(x  - x2  + 4xAs  + 4 xAs2) 


v _2 


+ g'  (0g"(.)  (x  + 2xAs) 


Again,  we  use  the  simplifying  assumptions  (see  Appendix  A. 2) 


3 
x 

4 
X 


- v2 
3 x 


Then,  (a, 13)  becomes: 

(*•14)  vx  = g'(.)2  I + V'(.)2(x2  + 2xAs2) 
+ 2g'(.)g"(.)  xAs 
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In  summary,  we  have: 


vm  * g(.)  + Jig"(.)  n 

v . 2 v 2 v 2 

\ “ g*(.)  n + >SSM(.)  nZ 

vx  - g(.)  + g’(.)As  + Jsg"(.)  (x  + Ai2) 

V = g’(.)2  X + Jsg"(.)2  (x2  + 2xAi2) 


+ 2g,(.)g"(.)  xAs 


These  equations  are  shown  as  Equations  (2.7)  through  (2.10), 


By  subtraction,  we  get: 


(vx  - vm)  = g'(.)As  + Jjg"(.)  (As  +As2) 


(vx  " V = g'(.)2As  + lsg 


;"(.)2£? 


v V V V 

As(x  + n)  + 2xAsz 


+ 2g?  (. )g"(.)  xAs 


These  equations  yield  the  correction  terms  for  the  effect  of  the 
Distribution  Approximation  and  are  shown  as  Equations  (2.11)  and 
(2.12). 
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Appendix  A,  2;  Simplifying  Assumptions 

In  Chapter  2,  we  assessed  correction  terms  for  the  mean  and 
variance  of  the  state  variable  and  calculated  the  resulting  correc- 
tion terms  for  the  mefin  and  variance  of  the  outcome  variable.  For 
several  linked  submodels,  we  saw  that  the  output  variable  of  one 
submodel  becomes  the  input  variable  of  the  next  and  that  we  can 
sequentially  calculate  the  correction  terms  for  each  intermediate 
variable. 

Clearly,  since  we  want  to  deal  only  with  the  mean  and  variance 
of  each  variable,  we  must  eliminate  all  terms  of  higher  order.  In 
developing  the  correction  terms  for  the  effect  of  the  Distribution 
Approximation  in  Appendix  A.l,  we  encountered  the  third  and  fourth 
central  moments  of  the  input  variable  and  eliminated  them  by  assuming 
the  following: 

(s  - s)3 
(s  - s)4 

If  we  know  only  the  mean  and  variance  of  a variable,  we  cannot 
reliably  estimate  its  third  central  momsnt,  which  is  a measure  of  the 
skewness  of  the  distribution  about  the  mean.  We  cannot  even  predict 
its  sign.  Therefore,  in  the  methodology,  I make  the  assumption  that 
the  third  central  moment  is  equal  to  zero,  which  is  true  if  the  dis- 
tribution is  symmetric. 


, v2 
3 s 
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On  the  other  hand,  we  can  estimate  the  fourth  central 
moment  if  we  know  the  variance-,  both  are  measures  of  the  dis- 
tribution's dispersion  about  the  mean.  In  fact,  it  turns  out 
that  the  fourth  central  moment  is  proportional,  or  nearly  so,  to 
the  square  of  the  variance  for  many  of  the  named  distributions. 
(See  Table  A.l)  The  ratio  of  the  fourth  central  moment  to  the 
variance  squared  ranges  from  1.8  for  the  uniform  distribution 
to  9 for  the  exponential,  but  many  of  the  distributions  have 
ratios  of  about  3.  Therefore,  I assume  for  the  purposes  of  the 
methodology  that  the  fourth  central  moment  is  equal  to  three  times 
the  variance  squared. 


Table  A.l 


Distribution 

(parameters) 


Ratio  of  the  fourth  central 
moment  to  the  variance  squared 


1.  Uniform  (a,b) 


a b 

A.  Normal  (y,a) 


9 

5 


for  any  a,b 


a 

2.  Triangular 

/T 

b 

(a.b.c) 

12 

5 

for  any  a,b,c 

a b 

3.  Trapezoidal 

c 

(a,b,c) 

12 

(1+8c+9c2+8c3+cA)  f 

h 

r 

|ch 

5 

234  any 

(l+8c+18c  +8c  +c  ) a,b 

for  any  y,o 


5.  Beta  (a,b) 


, (b+1)  /«  2b 

3 (b+2)  (b+3)  |_v  ' a(a-b) 


6. 

Beta  with  b - 2 a 

3 Iktil 

(b+3) 

7. 

Exponential  (a) 

9 

for 

any 

8. 

Laplace  (X) 

6 

for 

any 

9. 

Gamma  (a, 8) 

3 + — 

a 

for 

any  8 

10. 

Binomial  (n,p) 

3 6 4- 

1 

3 - — + 
n 

n2p(l-p) 

11. 

Poisson  (y) 

3 + ? 
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This  appendix  supplements  Section  2.3.  pp.  40-41.  We 
characterize  the  Dependence  Approximation  by  e : 

(b.l)  v ■ v + e 
a x 


Taking  the  mean: 

(b.2)  v •=  v + e 

a x 

Squaring  (b.l)  and  taking  the  mean: 


(b.  3) 

2 2 . 
v - v + e 

a x 

2 

1 + 2v  e 

X 

Squaring 

(b.2)  : 

(b.4) 

-2  - 2 
v - v + 

a x 

-2  - - 
e + 2v  e 

X 

Subtracting  (b.4)  from 

(b.3)  : 

(b.5) 

V V V 

v ■ v + e 

ax 

+ 2cov(v 

X 

So,  from  (b.2)  and  (b.5)  , we  have: 


v V v 

(va  " vx)  " e + 2cov(vx,r) 

These  equations  yield  the  correction  terms  for  the  effect  of  the 
Dependence  Approximation  and  are  shown  as  Equations  (2.14)  and  (2.15) 


Appendix  A. 4:  Calculation  of  e , e and  cov(v.e) 


This  appendix  supplements  Section  2.3,  p.  42.  We  can  calculate 
v 

e , e and  cov(vx»e)  from  the  assessed  conditional  mean  and 
variance  of  e , <e|d,s,4)  and  V(e|d,s,4)  as  follows: 


(c.l) 


(c.2) 


<e|d,4> 

j <e|d,s,4){s j4} 

J s 

<e2|d,4> 

j <e2 |d,s,4)  {s |4} 

•*  s 

( v<e|d,s,4)  + <e|d,s,4>2){s  1 4 } 

•*  s 

V<e|d,4> 


2 -2 
e - e 


( V<e|d,s,4>  + (e|dJs,4>2){s|4} 

* a 


<\e  |d,4> 


! [ g(d,s)  e {s,e|d,4} 

J6  Je 

j g(d,s){s|d,4}  e{e|d,s,4} 

s e 

g(d,s)  <e|d,s,4){s |4} 

* a 
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g(d,s)  <e|d,s,4){s|A}  - v^e 


Equations  (c.l)  , (c.2)  and  (c.3)  yield  the  required  calcula- 

tions and  are  shown  as  Equations  (2.16)  through  (2.18). 
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Appendix  A. 5;  Special  Caae  of  the  Dependence  Approximation 


Th^.s  appendix  supplements  Section  2.3,  pp.  44-45.  Suppose 
that  ^e  assess: 

<e |d»v  »4>  - A v 


Then: 


V<e|d,vx,4>  “ E vx 


<e|d,4>  - <e|d,vx,4>{vx|d,4} 

* 


A v {v  d,A} 

X X 


(d.l) 


'v 
A v 


e2  •»  <e2|d,A>  - <-i  |d,v  ,A>{v  |d,A> 


< V<e|d.".4>  + <e|d,v  ,4>2){v  |d,4} 

X XX 

V 

(B  vx2  + A2vx2){vjd,*} 


2 2 
* (B  + A ; v ^ 
x 


V 

e 


V<ejd,4> 


e - e 


*2n  2 .2-  2 
(B  + A ) v - Av 

X X 


(d.2) 


2 v - 2 

(B  + A ) v + B v 
x x 
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v e 
x 


(d.3) 


v 1 


V e{v  ,e|d,4} 
x x ' 


f vx{vx|d,4}je{e|d,vx,4} 


vx{vxld»4Keld*vv»A> 


A v iv  | d , A } 
X x' 


A v 


cov(vx,e) 


v e - v e 
X X 


. 2 - 2 

A v - A v 

X X 


V 

* A v 


Then,  the  correction  terms  for  the  effect  of  the  Dependence 
Approximation  are: 


(v  - v ) ■ e 

a x 


A v 


(v  - v ) ■ e + 2cov(v  ,e)  ■ (B  + 2A+A^)  v + B v 

ax  x xx 


These  equations  are  shown  as  Equations  (2.22)  and  (2.23). 
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Appendix  A.  6:  Multivariate  Case 

This  appendix  supplements  Section  2.5,  pp.  50-51.  Assume  that 
there  are  N state  variables.  We  use  the  following  matrix  notation: 

s_  , n and  21  are  N-dimensional  column  vectors 


w ■ Vg(d,s)  - , an  N-dim.  row  vector 

o * n 


W “ w w , an  N x N matrix 


i “ an  N-dim.  column  vector  of  all  l's 


L “ 1.  !_  , an  N x N matrix  of  all  l's 


($2g(d,s) 


'ij 


68^68^ 


s * n 


an  N x N matrix 


• l>  ] - 

:a  - Jcytiij  - SjXnj  - Sjlj  , 

;x  - [oyUt-IjXXj  -Sjjj  , 

- [Gljcov<n1,nJ)j 

- [OyCovCXj.Xjjj 

[GljV5j]  • 


an  N x N matrix 
an  N x N matrix 


G:cov(n) 
G:cov(jc) 
G:As  - 


, an  N x N matrix 
, an  N x N matrix 
an  N x N matrix 


Note  that,  for  a matrix  M , 1TM  1 is  the  sum  of  all  of  the 
elements  of  M. 
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f 


Also  note  that: 

- x - T 

(n  - n)  G (n  - n)  ■ 1 G:n  1^  , a scalar 

G:n  * G:cov(n) 

G:x  - G:cov(x)  + G:Aa_ 

We  characterize  the  Distribution  Approximation  by: 

As_  - x - n 

V V V 

As  ■ x - n 

We  expand  v - g(d,n)  about  n : 

m — — 

(e-D  = g(d,n)  + w(n  - n ) + % l1  G:n  1_ 

Taking  the  mean: 

(®*2)  = g(d,n)  + h G:cov(n)  1^ 

Squaring  (e.l)  and  taking  the  mean: 

vm2  = g(d,n)2  + lTW:nl  + UT  G:n  L G:n  1 

+ 2g(d,n)  w (n  - n)  + g(d,n)  1T  G:n  1 
+ w (n  - n)  lj  G:n  1^ 

(e*3)  vm2  = g(d,n)2  + W:cov(n)  1^  + h lj  G:n  L G:n  1 

+ g(d,n)  l1  G:cov(n)  1^  + w(n  - n)  G:n  1 
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Squaring  (e.2)  : 


_ 2 - 2 t 

(e,4)  vm  = g(d»n)  + h 1.  G:cov(n)  L G:cov(n)  JL  , 

+ g(d,n)  ij  G:cov(n)  1^ 

Subtracting  (e.4)  from  (e.3)  : 

(e,5)  \ ~ -T  W:cov(n)  I 

^ ^ I I^Gjn  l Gin  ■ G;cov(ii)  L G:cov(n)^  1 
+•  w (n  - n)  1T  G :n  1 

Now,  we  expand  vx  - g(d,x)  about  n : 

(e.6)  vx  a g(d,n)  + w (x  - n)  + 1 T G:x  1 

Taking  the  mean: 

(e.7)  vx  « g(d,n)  + vA£  + ^ 1T  |g:cov(x)  + G:Aij  1 

Squaring  (e.6)  and  taking  the  mean: 

vx2  ~ g(d,n)2  + JLT  W:x  1_  + \ 1T  G:x  L G:x  1 

+ 2g(d,n)  w (x  - n)  + g(d,n)  1T  G:x  1 
+ w (x  - n)  I1  G:x  1 

] 
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(e.8) 


g(d,n)2  + 1T  |w:cov(x)  + W:AsJ  1 


+ H ll  G:x  L G:x  1 

+ 2g(d,£)  wA£  + g(d,n)  1T  I G:  cov(x)  + G :t 


+ w (x  - n ) 1T  G:x  1 


G:cov(x)  + G:As  1 


Squaring  (e.7)  : 
(e.9) 


vx2  - g<d,n)2  + 1T  W:A£  1 


+ *s  1T  £g:cov(x)  + G:Asj  L |g:cov(x)  + G:AsJ 
+ 2g(d,n)  wAs  + g(d,n)  1T  j^G:Cov(x)  + G:As  j ] 

+ wAs_  l1  £g.*cov(x)  + G:Aij  1_ 

Subtracting  (e.9)  from  (e.8)  : 

(e.10)  vx  a 1T  W:cov(x)  1 + h,  1T  G:x  L G:x  1 

" * IT  |g:cov(x)  + G : Ai  j L £g:cov(x)  + G.-A^J  1 

1T  [gu 


+ w | (x  - £)  1T  gTx  1 - as 


cov(x)  + G:As_ 

Subtracting  (e.2)  from  (e.7)  and  (e.5)  from  (e.10)  , we 

obtain  the  estimated  correction  terms  for  the  effect  of  the  Distri- 
bution Approximation: 


(v 


vm)  ~ wAs  + H 1 1 J"g:cov(x)  - G:cov(n)  + G:A£ 


(e. 11) 


V V 

K - V 

x m 


(e.12) 


|w:cov(x)  - W:cov(n)^J  1^ 

+ H 1T  | G:x  L G:x  - G:n  L G:n 

- j^G:cov(x)  + G : A J L £g:cov(x)+G:AsJ 
+ G:cov(n)  L G:cov(n))  1^ 


+ w | (x  - n)  G:x  1_  - (n  - n)  1_T  G:n  1 
- As_  1^  |g:cov(x)  + G:  AsJ  l\ 


These  equations  are  not  useful  because  they  require  more  information 
than  we  have  about  the  two  distributions  on  b.  However,  when  the 
state  variables  are  mutually  independent,  these  equations  reduce  to 
a more  useful  form.  Letting  g^  and  g„  denote  the  partial  de- 

3 4 

rivatives  of  g and  letting  s and  s denote  the  third-  and 
fourth  central  moments,  we  have  for  mutually  independent  state 
variables : 


(e.13) 


l g1(d,n)As1  + \ ^g11(d,n)As1 
i i 

+ **£  i g^d.t^A^Aij 
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v _ 2 v 

(v  - v ) * £g,  (d,n)  As. 

x m 1 


<5  4 v 9 **  v 2 ^ — v - 2. 

+ H Ig11(d,n)2(xrx1  -n±+  n±  +4x^8^  4x1As1  ) 

i 

_2vv  vv  v 2,v  2, 

v “ ✓ . _ __  X.*  A r*  Xv  Ac  ' 


(e. 14) 


_ o v v vv  v - 2 v .-2* 

+ % I l g^W.n)  (xtx  -nin  -HXiAs^  +xjAs1  ) 
i J*i 

+ III  gi1(d.n)glk<d.n)(»i1Ai  + ^A^A^) 

i J WJ 


J J ^ 

+ J g1(d,n)g11(d,n)(x1  - h±  + 2x1As1) 

i 

+ J I g1(d,n)g1.(d,n)  2x±As 
i j*i 

To  eliminate  the  third-  and  fourth-order  terms,  we  apply  the 
simplifying  assumptions  (see  Appendix  A.  2)  : 


a 

s 

4 

s 


3 B2 


Then  (e.14)  becomes: 


v v 
(v  - v ) 
' x m 


l g1(d,n)2As1 


_ovv  vv  v _ 2 v - 2. 
+ *5  I I gyW.S)  ^xlx j ”nin j^^A® j + xjA8i  > 


(e.15) 


+ 1 I X8ij(d*S)8ik(d*S)  xiAijAik 

1 i 


+2 l l g, (d,n)g..(d,n)  x As 
i j 1 J J 


Equations  (e.13)  and  (e.15)  yield  the  estimated  correction  terms 
for  the  effect  of  the  Distribution  Approximation  for  independent  state 
variables  and  are  shown  as  Equations  (2.27)  and  (2.28). 
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